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Abstract 

The paper first discusses theoretically the off-resonance selective excitation 
method that is dependent on the atomic internal states and used to generate 
approximately a standard coherent state of harmonic oscillator. The coherent 
average method then is proposed to construct the state-selective trigger pulse. 
A state-selective trigger pulse can keep Gaussian shape unchanged but change 
in an internal-state-dependent form the center-of-mass position and/or momen- 
tum of an atomic Gaussian wave-packet motional state. A Gaussian wave-packet 
state is one of the simplest wave-packet states that can be easily manipulated 
and controlled in time and space. The paper also investigates how to manip- 
ulate in time and space an atomic Gaussian wave-packet motional state by a 
generalized quadratic potential field. A general quadratic Hamiltonian can af- 
fect not only the center-of-mass position and momentum but also the complex 
linewidth of a Gaussian wave-packet motional state while keep Gaussian shape 
of the motional state unchanged. It is shown that generally quadratic terms of a 
quadratic Hamiltonian can control directly the complex linewidth, while linear 
terms of a quadratic Hamiltonian can affect only the center-of-mass position 
and momentum of a Gaussian wave-packet motional state. 

1. Introduction 

In the paper [1] a particle picture has been used to describe intuitively how 
the halting-qubit atom evolves in the state-locking pulse field and how to con- 
struct a quantum control process to simulate the reversible and unitary halting 
protocol that is insensitive to its input state, although the quantum mechanical 
wave-packet picture has also used extensively in that paper. A particle picture 
is particularly intuitive to describe the decelerating and accelerating processes 
of a free atom and the elastic collision process for an atom bouncing off a hard 
potential wall. The correctness of a particle picture is based on the fact that 
a particle picture is very close to a wave-packet picture in quantum mechanics 
[2]. However, a quantitative and exact calculation in quantum mechanics for 
the time evolution process of a quantum system such as an atom does not use 
the classical particle picture. One must use wave-packet states or more gen- 
erally quantum states to calculate the time evolution process of the atom in 
quantum mechanics. The wave-packet states were used by Schrodinger, Dirac, 
and others to describe the quantum mechanical behavior of a particle date back 
to the early time of the wave mechanics. A very familiar example is that a free 



particle such as an atom may be described exactly by a Gaussian wave packet 
state in quantum mechanics [2]. Wave-packet states and especially Gaussian 
wave-packet states have been used frequently to describe the quantum collision 
and scattering processes in the atomic and molecular systems [2, 3, 4]. It can 
be seen in Refs. [5, 6, 7] that there is a more extensive application of the Gaus- 
sian wave-packet states to describe a variety of quantum dynamical processes 
of the atomic and molecular systems. Though quantum mechanically it is not 
limited to use the wave-packet states to describe the quantum control process 
of the reversible and unitary halting protocol and to investigate the mecha- 
nism of the state- locking pulse field [1], the quantum mechanical wave-packet 
states not only provide an intuitive picture for understanding the mechanisms 
of the reversible and unitary state-insensitive halting protocol and the state- 
locking pulse field but also simplify the quantitative and exact calculation for 
the quantum control process in the physical system of the halting-qubit atom. 
According to the reversible and unitary state-insensitive halting protocol [1] 
the wave-packet motional state of the halting-qubit atom should have a small 
spread, the wave-packet amplitude of the motional state decays quickly with 
the deviation from the wave-packet center and it is close to zero outside the 
effective spread. Obviously, such a wave-packet picture is close to a particle 
picture, which is just required by the reversible and unitary state-insensitive 
halting protocol. A Gaussian wave-packet state satisfies this requirement very 
miich. It can have a very narrow linewidth or wave-packet spread and its am- 
plitude at a position deviating from the center-of-mass position decays rapidly 
and exponentially with square of distance between the position and the center- 
of-mass position. A Gaussian wave-packet state is described completely by the 
three basic parameters: the center-of-mass position, the mean momentum, and 
the complex linewidth [2, 5, 6, 7, 8]. Then a Gaussian wave-packet state is 
simple and easy to be treated in theory and it is also easily manipulated and 
controlled in experiment. On the other hand, it is well known that the ground 
state of a harmonic oscillator is a Gaussian wave-packet state [2] , while at the 
initial time of the quantum control process the halting-qubit atom is prepared 
to be in the ground state of harmonic oscillator [1]. Thus, it is convenient 
and natural to choose the Gaussian wave-packet state to describe and calculate 
quantitatively the quantum control process. Due to the fact that a Gaussian 
wave-packet state has these advantages it should be better to keep Gaussian 
shape of the motional state of the halting-qubit atom unchanged in the whole 
quantum control process. 

Unitary manipulation in time and space in a quantum system plays a key 
role in implementing the reversible and unitary state-insensitive halting pro- 
tocol and realizing the efficient quantum search process [1, 12]. It is known 
that the quantum control process to realize the reversible and unitary state- 
insensitive halting protocol is a unitary evolution process in time and space [1]. 
One of the key components of the quantum control process is the time- and 
space-compressing processes which are realized by the unitary decelerating and 
accelerating processes of the halting-qubit atom. It has been shown [8] that 
the unitary decelerating and accelerating processes for a free particle moving 



2 



in space can be realized by the stimulated Raman adiabatic passage (STIRAP) 
method [9] . In the ideal or near ideal adiabatic condition the STIRAP pulse se- 
quence can transfer completely one Gaussian wave-packet motional state of the 
free atom to another [8] . It is known that the unitary propagator of a quadratic 
Hamiltonian can generally keep Gaussian shape unchanged for a Gaussian wave- 
packet state when the state is acted on by the unitary propagator [5, 6, 13, 14]. 
The Hamiltonian of the atom in the presence of the STIRAP pulse sequence 
is completely different from a conventional quadratic Hamiltonian. But it is 
surprising that the unitary propagator corresponding to the Hamiltonian of the 
atom in the presence of the STIRAP pulse sequence does not yet change Gaus- 
sian shape of an atomic Gaussian wave-packet motional state when the atom 
is irradiated by the STIRAP pulse sequence in the ideal or near ideal adiabatic 
condition. The STIRAP-based unitary decelerating and accelerating processes 
can manipulate the center-of-mass position and momentum of a Gaussian wave- 
packet state of a free atom in time and space, but it generally does not control 
the complex linewidth of the Gaussian wave-packet state because the imaginary 
part of the complex linewidth always increases linearly with the time period of 
the decelerating or accelerating process. One advantage to use the laser light 
fields such as the STIRAP pulse sequence to manipulate an atomic Gaussian 
wave-packet motional state is that the space-selective manipulation of the atom 
can be implemented easily. On the other hand, a general quadratic Hamiltonian 
can be used to manipulate not only the complex linewidth but also the center- 
of-mass position and momentum of a Gaussian wave-packet state. This will be 
investigated in detail in the paper. 

This paper is devoted to the construction of the state-selective trigger pulse 
and the unitary manipulation of a Gaussian wave-packet state in time and space 
by using a general quadratic Hamiltonian. A state-selective trigger pulse is a 
key component to realize both the reversible and unitary state-insensitive halt- 
ing protocol and the efficient quantum search process based on the unitary 
quantum dynamics in time and space [1, 12]. It is generally related to unitary 
manipulation and control in time and space of the atomic center-of-mass and 
internal motions as well as the coupling of the two motions. In the quantum 
control process the state-selective trigger pulse is used to transfer the ground 
motional state of the halting-qubit atom in the harmonic potential field to a 
standard coherent state [10, 11] with a higher motional energy. It is known that 
a standard coherent state of harmonic oscillator is a Gaussian wave-packet state 
[13]. There is a requirement that the coherent-state excitation process induced 
by the state-selective trigger pulse be dependent on the atomic internal state 
and keep Gaussian shape of the atomic motional state unchanged. One of the 
convenient methods to construct the state-selective trigger pulse could be that 
the Hamiltonian of the halting-qubit atom in the presence of state-selective trig- 
ger pulse is prepared to be a state-dependent quadratic Hamiltonian. A general 
quadratic Hamiltonian of an atomic system may be generated either by an exter- 
nal static electric or magnetic field or by the externally applied electromagnetic 
wave field. An external static electric or magnetic field usually does not induce 
the transitions of the atomic internal electronic or nuclear spin states. There- 
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fore, these uniform static electric and magnetic fields arc not suited to build 
up the internal-state-selective trigger pulse. Generally, the construction of the 
state-selective trigger pulse needs to use the external electromagnetic wave field. 

2. The constructions for the state-selective trigger pulse 

The state-selective trigger pulse is one of the key components to constriict the 
quantum control process to simulate the reversible and unitary state-insensitive 
halting protocol. It is used to control the halting operation in the reversible and 
unitary state-insensitive halting protocol. The properties for the state-selective 
trigger pulse have been described in detail in the paper [1]. There are two 
basic properties for the state-selective trigger pulse. One of which is that the 
trigger pulse is intcrnal-statc-sclectivc. Only when the halting-qubit atom is in 
the given internal state can the state-selective trigger pulse take a real action 
on the atom and excite the ground motional state of the atom to the coherent 
state. Another is that the state-selective trigger pulse can keep Gaussian shape 
unchanged for the atomic motional state. Under the state-selective trigger pulse 
the initial motional state of the atom could be converted into a standard coher- 
ent state of harmonic oscillator with a higher motional energy. In the quantum 
control process the initial motional state generally is prepared to be the ground 
state of the harmonic oscillator — the halting-qubit atom in the left-hand har- 
monic potential well [1] . It is well known that the ground state of a harmonic 
oscillator is a Gaussian wave-packet state [2] and so is a standard coherent state 
of the harmonic oscillator [13]. Therefore, it is clear that the state-selective trig- 
ger pulse should generate a unitary propagator for the halting-qubit atom such 
that the unitary propagator does not change the Gaussian wave-packet shape of 
the motional state of the halting-qubit atom. According to these properties one 
may construct explicitly such a state-selective trigger piilse for the reversible 
and unitary halting protocol in the physical system of the halting-qubit atom in 
the harmonic potential well. One could choose a sequence of laser light pulses as 
the state-selective trigger pulse in the atomic physical system. But there could 
also be other schemes to construct the state-selective trigger pulse. 

2.1. The non-internal-state-selective excitation process 

Before the internal-state-selective excitation process of an atom in a har- 
monic potential field is discussed in next subsections, this subsection is devoted 
to the investigation for the non-internal-state-selective excitation process of the 
atom. If there is not the internal-state-selective requirement, it is usually easy 
to excite a coherent state from the ground state of a harmonic oscillator. In 
principle one of the simplest methods to generate non-selectively a coherent 
state from the ground state of a harmonic oscillator is to apply directly an ex- 
ternal classical electric or magnetic field to the harmonic oscillator, so that an 
extra force generated by the external electric or magnetic field is exerted on 
the harmonic oscillator. Then it can prove [14, 15, 16, 17] that this extra force 
can drive the atom from the ground state to the coherent state with a higher 
motional energy. For example, for a trapped atomic ion in the harmonic poten- 
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tial well one may apply a spatially uniform driving electric field to the trapped 
atomic ion. The experiments to confirm the simple method have been carried 
out in the trapped ion systems [18, 19, 20]. Since the classical driving field does 
not excite any internal state of the trapped atom in the harmonic potential well 
and the trapped atom in any internal state can be excited from the ground 
motional state to the coherent state by the external driving field, this method is 
an intcrnal-statc-indcpendent excitation method to generate the coherent state. 
The coherent-state excitation process may be described simply by the Hamil- 
tonian of the atom in the harmonic potential well and in the presence of the 
external driving field, 

H{t) = Ho + H{r)+Hi{t), (1) 

where the Hamiltonian Hq describes the center-of-mass motion of the atom in 
the harmonic potential well and in the absence of the external driving field, 

here the oscillatory frequency lu = Lo{t) of the harmonic oscillator may be time- 
dependent or may not, the term H(r) is the internal Hamiltonian of the atom 
which describes the internal electronic (or nuclear spin) motion of the atom, 
and Hi{t) is the interaction between the atom in the harmonic potential well 
and the external driving field. One may not consider the internal Hamiltonian 
II{r) of the atom in the internal-state-independent excitation process, since the 
excitation process is not dependent on any internal state of the atom. One of 
the simplest interactions is Hx{t) = f{t)x [14, 15, 16, 17]. It indicates that there 
is the force exerted on the harmonic oscillator to generate the coherent state. 
This linear interaction Hi (t) may be generated by simply applying the driving 
electric or magnetic field to the harmonic potential well in a similar way to gen- 
erating the harmonic potential For example, if the atom is a charged 
atomic ion with charge q, the interaction may be written as Hi(t) = —qxE{t) for 
the spatially uniform and time-dependent driving electric field E{t) and hence 
the forced function f{t) is given by f{t) = —qE{t). The atom is usually in some 
hypcrfine ground electronic (or nuclear spin) state before the external driving 
electric field is applied. However, when the atom is irradiated by the suitable 
external electromagnetic field, it may jump to an atomic excited state or an- 
other hyperfine ground state. These internal electronic (or nuclear spin) state 
transitions of the atom have their own transition frequencies. Suppose that the 
minimum transition frequency among these internal-state transition frequencies 
is much greater than the oscillatory frequency of the external driving electric 
field E{t). Then it is impossible for the external driving electric field to in- 
duce the atom to make a transition between the atomic internal states. During 
the external driving field the atom stays in its initial internal state. Thus, the 
coherent-state excitation process of the atom under the external driving electric 
field is independent of the atomic internal states and hence it is internal-state- 
independent. On the other hand, while the coherent-state excitation process is 
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independent of any atomic internal state, the driving electric field E{t) must be 
able to induce effectively the transitions between the vibrational energy levels of 
the atom in the harmonic potential well. The vibrational energy levels generally 
have a much smaller energy-level space than the atomic internal energy levels. 
If now the oscillatory frequency of the external electric field is set to be near 
the oscillatory frequency of the harmonic oscillator — the halting-qubit atom in 
the harmonic potential well, then the external driving electric field may induce 
the atom to make a transition between the vibrational energy levels, while the 
atomic internal states are not affected by the driving electric field. Since any in- 
ternal state of the atom remains unchanged during the coherent-state excitation 
process by the external driving electric field, the atomic internal Hamiltonian 
H{r) may be omitted from the total Hamiltonian H{t) of Eq. (1). Then in this 
case the Hamiltonian to describe the time evolution process of the atom in the 
harmonic potential well and in the presence of the external driving electric field 
is reduced to the form 

H{t) = 1^ + ^mLu^x^ + f{t)x. (3) 

This is the Hamiltonian of a time-dependent forced harmonic oscillator [13, 14, 
15, 16, 17, 21]. If now the initial motional state of the atom is the ground state 
of the harmonic oscillator, then the unitary propagator of the Hamiltonian H{t) 
(3) will convert the ground state into a coherent state [17], indicating that the 
external driving electric field can transfer the ground state to a coherent state. 
Though the construction for the state-selective trigger pulse does not use the 
state-independent excitation method as mentioned above, this method may be 
useful when one manipulates unitarily a Gaussian wave-packet state in time 
and space by a unitary propagator of a quadratic Hamiltonian, as can be seen 
in section 4. This is because the Hamiltonian (3) is really the specific form 
of a general quadratic Hamiltonian. The internal-state-independent excitation 
method is also useful in the quantum control process [1]. When the halting- 
qubit atom returns back to the left-hand harmonic potential well in the quantum 
control process, its motional state needs to be transferred back to the original 
ground motional state. Then the state-independent excitation method may help 
the quantum control process to realize such a state transfer. 

For a neutral atom system both the harmonic potential field and the external 
force field may be generated by applying an external driving electric or magnetic 
field. It is known that an atom can have an induced electric dipole moment in the 
presence of an external electric field [2]. Then a potential energy V{x^t) of the 
atom can be generated due to the induced electric dipole moment in the external 
electric field E{x, t), which may be expressed as V{x, t) = ^a\E(x, where a 
is the atomic polarizability [2]. And hence the interaction Hi{t) in Eq. (1) could 
be obtained from the potential energy V{x,t). Obviously, this atomic potential 
energy is internal-state-independent. Thus, the time-dependent Hamiltonian 
H(t) of Eq. (3) for the neutral atom in a harmonic potential well may also 
be generated by applying a suitable time- and space-dependent electric field 
E{x, t) to the atom. For example, if the external electric field E{x, t) oc ax + b, 
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then the potential energy V{x,t) contains the hnear interaction Hi(t) oc x. 
On the other hand, a nuclear spin or a neutral atom which has an intrinsic 
magnetic dipole moment may generate a state-dependent force in an inhomo- 
gcneous magnetic field [50]. A nuclear spin or atom in an external magnetic 
field can generate the Zeeman effect [2, 50] and the spin energy level is given 
by Em = —jB{x,t)fim, where B{x,t) is the external magnetic field strength, 
m, the spin magnetic quantum number, and 7 the gyromagnctic ratio of the 
spin. Thus, the interaction between the spin and the external magnetic field is 
given by Hi{t) = —^fiB{x,t)m. When the external magnetic field is not uni- 
form, for example, B(x,t) oc x, the spin will be acted on by an external force 
generated by the magnetic field B{x,t) if the spin magnetic quantum number 
m ^ 0. Since the quantum number m marks the spin quantum state, the spin 
state with m = is not acted on by the external magnetic field, while all those 
spin states with m ^ undergo the external magnetic field. Therefore, the spin 
interaction Hi{t) may be time- and state-dependent. Then an inhomogeneous 
external magnetic field could act as a state-selective trigger pulse if the halting- 
qubit atom is chosen as a nuclear spin (its spin quantum number is even) or 
a neutral atom with an intrinsical magnetic dipole moment and the external 
magnetic field is designed suitably. 

2.2. The off-resonance selective excitation method 

A general method to construct the state-selective trigger pulse is involved 
in manipulating and controlling in time and space the center-of-mass motion 
and the internal electronic (or spin) motion of an atom as well as the coiipling 
between the center-of-mass and the internal motion. Since the coherent-state 
selective excitation process is dependent upon the specific internal state of the 
atom during the state-selective trigger pulse, the atomic ccnter-of-mass motion 
must be coupled with the internal electronic (or spin) motion of the atom in 
the excitation process. Coherently manipulating the center-of-mass and the 
internal motion of the atom as well as the interaction between the two mo- 
tions is the fundament for realizing the atomic laser light cooling and trapping 
in an atomic ensemble [22], implementing the quantum computation [23], and 
preparing and transferring various quantum states in an atom-ionic system in a 
harmonic potential field [24, 25, 26, 27]. The state-selective trigger pulse is also 
closely related to the coupling between the atomic center-of-mass and internal 
motions. Laser light is a general technique to realize the coupling between the 
internal and the center-of-mass motion of an atom. In particular, as one of the 
most useful double-photon excitation methods the stimulated Raman adiabatic 
passage (STIRAP) method [9] has been extensively used to coherently manip- 
ulate the center-of-mass and the internal motion of an atom and also used to 
create and control the coupling between the atomic center-of-mass and internal 
motions in an atomic ensemble in quantum interference experiments [28, 29]. 
The largest advantage of the STIRAP method is that the STIRAP method can 
achieve a complete state-transfer efficiency in theory and is tolerant to the ex- 
perimental imperfection and can avoid the spontaneous emission generated by 
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atomic excited states. For example, a high internal-state transfer efficiency is 
achieved experimentally by the STIRAP method in trapped ions [30]. On the 
other hand, the conventional Raman double-photon laser light techniques also 
have been used to prepare the ground motional state and various quantum coher- 
ent states and to transfer one atomic internal state to another in an atomic ion 
system in a harmonic potential field [18, 27, 31, 32, 33]. The STIRAP method 
usually uses a pair of Raman adiabatic laser light beams to couple the atomic 
center-of-mass motional state and internal electronic states (or spin polarization 
states) so as to realize the interaction between the atomic center-of-mass and 
internal motions. In order to generate effectively and internal-state-sclcctively 
the coupling between the two atomic motions the laser light electromagnetic 
field should have an oscillatory frequency close to the resonance frequency of a 
given pair of the internal electronic (or spin) states of the atom and far from 
those resonance frequencies of any atomic internal-state transitions other than 
this given resonance frequency. With this frequency setting the electromagnetic 
field can induce only the transition between the given pair of the atomic inter- 
nal states and does not induce any other atomic internal-state transitions when 
the amplitude of the electromagnetic field is not very large. Thus, this is an 
internal-state-sclcctive excitation process. This state-selective excitation pro- 
cess may be described by the unitary quantum dynamics in theory. Since this 
is an internal-state-selective excitation process, one must consider the atom to 
be a multi-level physical system. The simplest case is that the atom is a three- 
level physical system and two of the three internal energy levels are irradiated 
selectively by the external electromagnetic field. Then in this case the atomic 
internal Hamiltonian H{r) may be simply written as 

H{r) = Eo\Mr)){Mr)\ + E^\^,{r)){Mr)\ + E2\^2{r)){Mr)\ (4) 

where Ek and \ipi,{r)) {k = 0, 1, 2) are the k—th eigenvalue and eigenstate of the 
internal Hamiltonian H{r), respectively, that is, H{r)\ipf.{r)) = Ek\ipki''')) ■ For 
example, the eigenstates |V'o(^)) ^^"^ IV'i('')) be taken as the two hyperfine 
ground electronic states \go) and [^i) of an atom, respectively, while \tp2i''~)) 
taken as the atomic excited state |e). Here suppose that the external electro- 
magnetic field is applied selectively to the two atomic internal energy levels j^o) 
and ]e), while the internal state [^i) is not affected by the electromagnetic field. 
On the other hand, the semiclassical theory of the electromagnetic radiation 
and the electric dipole approximation are still suited to describe the state- 
selective excitation process in the physical system of the atom plus the elec- 
tromagnetic field [2, 34]. In the electric dipole approximation the interaction 
between the atom and an externally applied electromagnetic field can be ex- 
pressed as Hi{t) = —D . E(a;, t), where D is the electric dipole operator of the 
atom, E{x, t) the time-dependent electric field of the externally applied electro- 
magnetic field, and the coordinate x the center-of-mass position of the atom. 
The selective excitation method may use either the single- or double-frequency 
(or double-photon) or even multi-photon excitation method. For the single- 
and double-frequency selective excitation processes the external electromagnetic 
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fields are respectively written as 

E{x, t) = ^ELo{t) exp[i{kLo-x - ULot)] + C.C. (5) 

and ^ 

^{^^ t) = 2^Lo{t) ex.p[i{kLo-x - LOLot)] 

+ ^ELiit) exp[i{kLi.x - ULit)] + C.C, (6) 

where ELk{t), kLki and uilh {k = 0, 1) are the complex amplitude, wavevector, 
and oscillatory frequency of the laser light beam, respectively, and C.C. stands 
for the complex (or Hermite) conjugate term. In the single-frequency selective 
excitation process only one laser light beam {ELo{t), fcio, wlo) is selectively ap- 
plied to the two internal energy levels \go) and |e), while in the double-frequency 
selective excitation process a pair of the laser light beams {ELo{t), kiOj^^Lo) and 
/cli, wli) are selectively applied to the two internal states \go) and |e) 
simultaneously. Note that the atomic internal state \gi) is not affected by any 
laser light beam in the single- and double-frequency excitation processes. There 
is a slight difference between the double-frequency selective excitation method 
here and the conventional STIRAP experiment. In the conventional STIRAP 
experiment a pair of the Raman laser light beams usually are applied to two dif- 
ferent transitions linking the two different ground internal states \go) and \gi) to 
the same excited state |e), respectively. Now the electric dipole interaction Hi{t) 
between the atom and the electric field E{x, t) of the external electromagnetic 
field is given explicitly by [24, 34, 35] 

ifl(i) = n{^L0{t)I+ + nio{t)I-)cOs{kLoX - UJLot) 

= Hflo{t){I~^ exp[i{kLox - uiLot - <^o(*))] 

+1' cxp[-i{kLox - ujLot - fo(t))]} (7) 

for the single-frequency selective excitation process, where the second equal- 
ity is obtained in the rotating wave approximation, and in the rotating wave 
approximation for the double-frequency selective excitation process, 

Hi{t) = Mlo{t){I~^ exp[i{kox — ujot — <^o (*))] 

+/" exp[-i(fcoa; - uot - (Po{t))]} 
+fiD,i{t){I~^ exp[i{kix — uit — <^i(i))] 

+I~ exp[—i{kix — u)it — (fii{t))]} (8) 

where the two laser light beams may be either counterpropagating (/co and fci 
have the opposite signs) or copropagating {ko and ki have the same signs), the 
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amplitude fluit) = ri;(t) exp[— iy>;(t)] {I = 0,1), and the atomic internal-state 
operators are defined by 

21, = (/I - /O), /° = \go){go\, = |e)(e|, /+ = \e){go\, /" = \go){e\. 

These interactions (7) and (8) are similar to those in the Jaynes-Cummings 
model [35] of an atom plus electromagnetic field system. Now the total Hamil- 
tonian (1) can be given explicitly if the internal Hamiltonian H{r) (4) and the 
electric dipole interaction Hi{t) of Eq. (7) or (8) are inserted into Eq. (1). It 
should be pointed out that it is possible to apply an extra laser light beam for 
each laser light beam above to compensate the rotating-wave approximation. If 
each laser light beam in the selective excitation process above is replaced with 
a pair of the laser light beams with the orthogonal electric field vectors and the 
suitable phases, one can eliminate the rotating-wave approximation. This means 
that the Hamiltonian (8) may be constructed exactly. 

It is known that the quantum behavior of an atom in the harmonic potential 
field and in the presence of the external electromagnetic field may be described 
by the complete set of the product basis states {\^nk{x,r))}, 

|*nfe(ar,r)) = |^„(a;))|Vfe(r)), (9) 

where [ijjni'-^)) is an eigenstate of the Hamiltonian Hq (2) of the harmonic oscil- 
lator, which is used to describe the atomic center-of-mass motion, while \tpk{r)) 
is an eigenstate of the internal Hamiltonian H{r), as can be seen in Eq. (4), 
which is used to describe the internal electronic (or spin) motion of the atom. 
Since the external electromagnetic field is applied to only the two internal states 
\go} and |e), the time evolution process of the atom in the harmonic potential 
well and in the presence of the external electromagnetic field is described by the 
unitary propagator 

U{t) = Texp{-1 dt'H{t')} = eM^i^\g,){g,\}UUt) (10) 

where the total Hamiltonian H{t) is given by Eq. (1) and the unitary propagator 
UL{t) is defined as 

UL{t) = TeM-^J^ dt'HLit')} (10a) 
with the Hamiltonian Hl (t) given by 

HL{t) = ^ + Imio^x^ + EoI° + E2l^ + Hi{t). (11) 

Though the Hamiltonian iJ^ [t) is involved in only the two-state subspace span 
by the internal states jgo) and |e), it contains the term iJo (2) and conse- 
quently the unitary propagator UL{t) still can affect the atomic product state 
\^{x,t))\gi), where \^{x,t)) is an atomic center-of-mass motional state. For ex- 
ample, the unitary propagator U l {t) may have a significant effect on the product 
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state t)) l^i) if the motional state \'^{x,t)) is a superposition state. Ac- 

tually, the Hamiltonian HL{t) could not be considered as the Hamiltonian of a 
two-level system consisting of the internal states l^o) and |e). However, the uni- 
tary propagator UL{t) can generate only a global phase factor for the product 
state |^'(x,t))|5i) if the motional state |4'(x,t)) is an energy eigenstate \^^{x)) 
of the harmonic oscillator. Now suppose that at the initial time the halting- 
qubit atom is in the product state \ijjQ{x))\gi) , where \iPq{x)) is the ground 
motional state of the harmonic oscillator. Obviously, the atom is still in the 
product state |V'o(a^))|g'i) if it is acted on by the unitary propagator U{t) of Eq. 
(10). However, after the atom is transferred to another internal state |go) or 
|e) from the initial internal state l^i), it could be excited to a motional state 
that has a higher motional energy than the ground motional state by the uni- 
tary propagator U{t). This means in the case that the imitary propagator U{t) 
is really state-selective, and the excitation process may be described directly 
by the Hamiltonian HL{t) of Eq. (11). Now one wants to design the external 
electromagnetic field E{x, t) of Eq. (5) or (6) such that the atomic ground mo- 
tional state is transferred to a standard coherent state of harmonic oscillator by 
the unitary propagator U{t) of Eq. (10) or Uhit) of Eq. (10a). The conven- 
tional Raman laser light beams [18, 19, 31, 32, 33] are often used to generate 
selectively the coherent state. They may be understood intuitively below. At 
the starting time of the excitation process the atom is in the ground motional 
state \i1^q{x)) and the internal state |go)- Then the first Raman laser light beam 
excites the atom from the internal state \go) to the excited state |e), while the 
second Raman laser light beam induces the atom from the excited state |e) to 
jump back to the original internal state jgo), and at the same time the initial 
ground motional state |^Q(a;)) is changed to a coherent state during the excita- 
tion process. Although the atomic internal state |g(o) is not changed after the 
excitation process, the atomic motional state is changed from the initial ground 
state \iPq{x)) to the coherent state. Obviously, this coherent state is indirectly 
generated through the internal-state transfer pathway jfifo) ~* |e) — > l^o) by the 
two Raman laser light beams. Therefore, this is an internal-state-dependent 
excitation process. 

It could be convenient to investigate the internal-state-selective excitation 

process based on the Raman laser light beams in the Heisenberg picture which 
is often used in the laser spectroscopy [34] . The excitation process is involved in 
only the two-state subspace span by the internal states [g^o) and |e) and is gov- 
erned by the Hamiltonian HL{t) of Eq. (11). Then all those operators appearing 
in the Hamiltonian HL{t) of Eq. (11) are defined as A{t) = UL{t)^AUL{t) in 
the Heisenberg picture, where the operator A may be 7°, 7^, 7+, 7~, x, and 
p. The dynamical equations for these operators in the Heisenberg picture are 
given by 

^hJ^A{t) = [A{t),HL^t)]. (12) 

Below consider a general double-frequency excitation method including those 
using the conventional Raman laser light beams. By inserting the Hamiltonian 
7fi(f) of Eq. (11) with the electric dipole interaction Hi{t) of Eq. (8) into the 
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Heisenberg equations (12) one obtains [34] 

= -tuvj^{t) + 2mo{t)I,{t) exp[-ifcoa;(i)] exp[i(wot + <^o)] 

+2nni{t)h{t) exp[-ifcia;(i)] exp[i(a;it + y>i)], (12a) 
ih^I~{t) = tiujal~{t) — 2Mlo{t)Iz{t) e-xp[ikox{t)] exp[— i(a;oi + <^o)] 

-2W.i{t)Iz{t) exp[ifcia;(t)] exp[-z(a;it + y-J], (12b) 
iH^hit) = Mlo{t){I~^{t) exp[i{kox{t) - uot - <^o)] 

-I~{t) exp[-i(fcox(t) - u)at - lpq)]} 
+nni{t){I+ {t) exp[i(fcia;(t) - ujit - </Ji)] 



and 



-I~{t) exp[-i(fcia;(i) - uj-^t - v?i)]}, (12c) 

^p{t) = -TmJ^x{t) 

—ihkQQ.o{t){I'^ {t) exp[j(fcoa;(t) — oj^t — <^o)] 

-I~ {t) eyc£>[-i{kox{t) - ujot - (^q)]} 
—ifikiQ.i{t){I^ {t) exp[j(fcia;(t) — uJit — <^i)] 

-I~{t) exp[-j(fcia;(i) - uJit - (12d) 

j^x{t)=p{t)/m, (12e) 

where EqI^ + i?2-f^ = ao-E^ + fii^alz, and _E is the 2x2 unit operator, = 
{E2 + Eo)/2, fvjJa = {E2 — Eq), and oja is the resonance frequency of the two 
atomic internal energy levels \go) and |e). The former three equations mainly 
describe the atomic internal motion and the coupling between the internal and 
center-of-mass motions of the atom, while the rest two equations mainly describe 
the center-of-mass motion and the coupling of the two motions. The Heisenberg 
equation set of Eqs. (12) describes completely the time evolution process of 
the atom, which is involved in the atomic internal and center-of-mass motions 
as well as the coupling of the two motions. In the absence of the external 
electromagnetic field the Heisenberg equations (12) have a simple solution [34]: 

I^{t) = 7±(0) e^p[±iujat], hit) = W), 

a{t)~^ = a(O)"'" exp(iwt), a{t) = a(0) exp(— iwf), 
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where the creation (a+) and annihilation (a) operators defined through 



""^^^^ V2^^"^*^''+"^*^^' Pit)=^^lf>^rn{a{t)+ -a{t)). (13) 

This simple solution is called the uncoupling solution to the Heisenberg equa- 
tions (12). 

Now the off-resonance excitation method is introduced below. The off- 

rcsonancc excitation means that the two atomic internal states |go) and |e) are 
irradiated by weak and off-resonance laser light beams. This also means that 
for the off-resonance excitation using the Raman adiabatic laser light beams 
the Rabi frequencies f7o(t) and ili{t) of the two Raman adiabatic laser light 
beams are slowly varying and much less than the frequency offsets {cOa — t^o) 
and {u)a — wi). Obviously, in the off-resonance excitation the solution to the 
Heisenberg equations (12) should be close to the uncoupling solution and hence 
the operators {I'^it) exp[^iu!at]} should be close to the operator /""^(O). There- 
fore, the time derivatives of the operators {I^{t) exp[^icOat\} are close to zero. 
Now one may make a unitary transformation [33] : 

i^{t) = exp{-icOaIzit)t)I^{t) exp{iijJz{t)t) = I^it) exp[TiWai]. (14) 

Then the first two Heisenberg equations (12a) and (12b) are reduced respectively 
to the forms 

^/+(t) = -2ifloit)Iz{t) exp[-ikox{t)] exp[-i{LOa - wo)* + iVo] 

—2iQi{t)Iz{t) exp[—ikix{t)] exp[— ^(wa — u}i)t + iifi)] (15a) 

and ^ 

— 7~(t) = 2iQoit)Iz{t) exp[ikox{t)] exp[i(a;o - ujo)t - iipo] 

+2iD,i{t)Iz{t) exp[ikix{t)] exp[i(a;a — uji)t — i(pi\. (15b) 

Obviously, the time derivatives of the operators {/^(t)} are close to zero when 
the Rabi frequencies f2o(i) and ili{t) are close to zero. If the frequency offsets 
{\{uja — t^;)!} much greater than the Rabi frequencies {ri;(t)} {I = 0, 1), the 
atomic motional velocity p{t)/m, is small in the harmonic potential well, and 
the time varying of the Rabi frequencies is slow (the Raman adiabatic laser 
light beams satisfy the condition), then the operators {I^{t)} may be obtained 
approximately by integrating the two equations (15), respectively, because the 
oscillatory terms exp[±i(a;; — ^^a)A [l ~ 0, 1) will generate a dominating con- 
tribution to the two integrals. Now by integrating first the two equations (15), 
then making integration by part, and then by using Eq. (14) one obtains the 
two operators: 

I+{t) - exp(iwai)/+(0) = '] ^'-''I'^'^'o) no{t)Iz{t) exp[-ikox{t)] exp[wot] 
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+ ^^^^^^ni{t)I,{t) exp[-ikix{t)] exp[i(Jit] + E{I+), (16a) 
r{t) - exp(-it^ai)/"(0) = '^f^^''~^'^{' no{t)h{t) exp[ikoxit)] exp[-iwoi] 

+ ^^^P(~''^i ) (t) (^) e^p[ikixit)] e^p[-iujit] + j;( J- ) (16b) 

(Wo — Wi) 

For simplicity, here suppose that the phases lpq and are time-independent and 
at the initial time both the amplitudes of the two Raman adiabatic laser light 
beams are zero, that is, ilo{0) = rii(O) = 0. Since these terms {exp(±iW(jt)/='=(0)} 
are of the uncoupling solution, the right-hand sides of the equations (16) mea- 
sure the deviation of the real solution from the uncoupling solution. The terms 
£'(/+) and E[I~) in Eq. (16) are error terms which are higher-order terms of 
the frequency offsets {{oja — toi)~'^ (/ = 0, 1), (wa — Wi)~^(wa — w;)""^, etc.). The 
error terms E(I^) and E{I^) can be neglected in the first-order approxima- 
tion. That is, only the first-order terms which are proportional to the inverse 
frequency offsets (jwo — Z = 0, 1) are retained in the operators {/^(t)} in 

the first-order approximation. Now inserting the operators {I^{t)} (16a) and 
(16b) of the first-order approximation into the Heisenberg equation (12c) one 
obtains, up to the first-order approximation, 

m - 4(0) = _Ml^E^ti^i+iQ) exp[*fcoa;(t)] exp[*(c.„ - u,,)i\ 

[Wa - Wo) 

-^^^^^^^^fc^/+(0)exp[ifcia;(t)]exp[i(cj„ -cji)f] +aC. (16c) 

(to a — Wi) 

where C.C. stands for the hermite conjugate of the first two terms on the right- 
hand side of Eq. (16c). The two Heisenberg equations (12d) and (12e) then can 
be solved with the help of the first-order approximation operators {I^{t)} (16a) 
and (16b) and Izit) (16c) and one obtains, in the first-order approximation. 



a{t)+ - exp(ia;t)a(0)+ = i/z(0)W — [ dt'{neff{t') 

X exp[-ia;(t' - t)] sin[Aa;t' -|- Atp]}, (16d) 



a{t) - exp{-'iujt)a{0) = -i4(0)W / dt'{neff{t') 



2H 







X exp[iui{t' - t)] sinfAwi' -|- A(p\}, (16e) 

where the Lamb-Dicke limit ||AA:a;(i)|| << 1 and the condition muP'/fi » 
|Afcf2e//(i)| have been used; AA: = ko — ki, Alo = loq — loi, Aip = (fg — (fi, and 
the parameter f2e//(i) is given by 



{(Ja-OJl) (Wa-Wo) 
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The condition muP >> |fiAfcf7e//(t)| means that the effect of the Raman laser 
Hght beams on the harmonic potential field (ma;^x^/2) can be neglected. Now 
by inserting the first-order approximation operators {/^(t)} of Eq. (16a) and 
(16b) into the electric dipolc interaction Hi{t) (8) and then inserting Eq. (8) 
into Eq. (11) the Hamiltonian HL{t) of Eq. (11) can be written as, in the 
first-order approximation, 

HlH) = aoE + tko[a+{t)a{t) + ^] + hna{t)I,{t) 

-|-2/iAfc[Oeo(i) + 0.^i{t)\lz{t)x{t) sinfAwf -|- A(p] 
-|-fif2o(t){exp(— i(/?Q)7"'"(0) exp[ikox{t)\ exp{i{wa — i^o)t] 

-|-exp(i(^o) exp[-jfcoa;(t)]7~(0)exp[-i(a;o - uJo)t]} 
-|-fif2i(t){exp(— iy>i)7"'"(0) exp[ikix{t)] exp[i(a;a — wi)t] 

+ cyiY>{i(pi) cx:p[—ikix(t)]I~ (0) exp[— i(wa — '^-'i)*]} (17) 

where the expansion exp[±zAA;a:(i)] = 1 ± il^kxit) + ... in the Lamb-Dicke limit 
has been used, the parameter Q.a{t) is given by 

^ait) =Ua + , + , ' . + 2[n,o{t) + n,,{t)] cos(Aa;t + A^), 

[UJa - Wo) {U)a - Wi) 

and the effective Rabi frequency i^ei{t) = 2Qo{t)^i{t) / {cJa — i^i) {I = 0, 1) which 
was also obtained in Rcfs. [36]. ft is known that the Lamb-Dicke parameter 77 
is defined through kx{t) = ri{a{t)'^ + a{t)). The Lamb-Dicke parameter r/; = 
y/{h'^kf)/{2mtkj) with the wave vector value ki = coi/c {I = 0, 1) may be large 
if the frequencies luq and uji of the two Raman laser light beams are near the 
transition frequency uja (~ 10^^) between the atomic internal states \go) and |e). 
However, the wave vector difference |Afc| = |fco — fci| may be very small and its 
value can be controlled in experiment. For example, |Afc| = |Aa;|/c ~ lu/c << 
oJa/c, since the oscillatory frequency of the harmonic oscillator ^ 10^ is much 
less than the atomic internal-state transition frequency uJa ~ 10^^. In this case 
the Lamb-Dicke parameter may be very small: rj = ^/ (Ak)"^ / {2mfkj) << 1, 
and hence ||Afca;(t)|| << 1. 

It seems that one may directly obtain the first-order effective Hamiltonian 
from Eq. (17) by omitting those terms of the operators {/^(O)} on the right- 
hand side of Eq. (17), but the correct treatment is that one should use the first- 
order approximation solution of the Heisenberg equations (12) which contains 
the operators /^(t), a{t)^, and a{t) of the equations (16c), (16d), and (16e) 
to reduce the Hamiltonian (17) to the first-order effective Hamiltonian. The 
Hamiltonian (17) may be divided into two parts, the first part does not contain 
the operators {/^(O)}, while the second contains. Note that the effective Rabi 
frequency flait) in the Hamiltonian (17) is close to uJa when the Rabi frequency 
fll{t) {I = 0, 1) is much smaller than \{uja — Wi)|. Then the part that contains 
the operators {/^(O)} is nonsecular, since its components generally contain 
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oscillatory phase factors cxp[±i{uja — ^i)t] {I = 0, 1) instead of exp[zLiuJit] (see 
Eq. (17)) and in the rotating frame defined by the Hamiltonian toilz or LOalz they 
still contain largely oscillatory phase factors exp[±iuat], exp[±^a;^^], etc., instead 
of the slowly varying phase factors cxp[±iALut\. Thus, the contribution of the 
nonsecular part to the off-resonance excitation is usually small. If one chooses 
the suitable parameters of the Raman laser light beams such that the absolute 
amplitudes in front of the oscillatory phase factors exp[±i(wa — LUi)t\ {I = 0, 1) 
in the nonsecular part are much smaller than | (wq — ) | , then the contribution 
of the nonsecular part to the off-resonance excitation may be neglected. For 
convenience to discuss below, here the parameters of the Raman adiabatic laser 
light beams are chosen suitably such that the first part is much more important 
than the second and hence the nonsecular part that contains the operators 
{/^(O)} may be neglected. Then in this case the Hamiltonian (17) is reduced 
approximately to the simple form 



where the unity operator term is omitted and the parameter fi{t) is given by 



In the parameter f3{t) the integral containing the parameter £l^ff{t) usually is 
not considered [36]. If the two frequencies offsets {ua — wj} {I = 0, 1) have 
the opposite sign to each other such that the parameter |f2e//(i)| is large, then 
the dominating contribution to the parameter /3(i) may come from the integral 
containing the parameter Qeffit). The first-order effective Hamiltonian (18) can 
be used to excite a coherent state of harmonic oscillator. The cross operator 
[P{t)a~^ + P{t)* a]Iz in the Hamiltonian HL{t) (18) is responsible for the internal- 
state-selective excitation of a coherent state. This will be interpreted below. In 
the interaction representation defined by the harmonic oscillator Hamiltonian 
Hq = lvjj{a^a + 1/2) the Hamiltonian (18) is transformed to [11] 



where the complex number z{t) = (3{t) exp[iu!t — i7r/2]. The propagator of the 
Hamiltonian HLi{t) in the interaction representation is 




(18) 




+Ak\ [Oeo(t) + f^ei(t)] sm[Aujt + Aip] exp{iujt). 




HLi{t) = ma{t)h + in[z{t)a+ - z{t)*a]h 




where the unitary displacement operator D{ii{t)) is defined as [11] 
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Note that here of the operator D{ii{t)) is really an operator and may be 
expressed /u(t) = iJLi{t)E + 2IzH2{i) pl], where {I = 1,2) is a complex 

number, E = (|e)(e| + |s'o)(S'o|), and 4 = (|e)(e| - \go){go\)/^- It is known 
that the imitary displacement operator D{z{t)) {z{t) is a complex number) can 
excite a standard coherent state from the ground state of a harmonic oscillator 
[11]. Since the operator ^{t) acts on only the two atomic internal states |e) and 
l^o), the imitary displacement operator D{ji{t)) can excite the coherent state of 
the atom in the harmonic potential well only when the atom is in the internal- 
state subspace span by the two internal states |e) and \go)- This indicates that 
the imitary propagator ULi{t) can excite the coherent state in an internal-state 
dependent form and hence is internal-state-selective. In order to see more clearly 
how the first-order effective Hamiltonian (18) can be used to excite selectively 
the coherent state one may assume that the parameter /3(f) is real, which could 
be achieved by choosing suitably the parameters of the Raman laser light beams. 
Then the Hamiltonian (18) is written as, with the help of Eqs. (13), 

HL{t) = ^ + Imu^x^ + Mla{t)h + f{t)hx (19) 

where the force function f{t) = \/2mfibjf3{t). The last term on the right-hand 
side of Eq. (19), which is proportional to the product operator J^x, is responsible 
for the internal-state-selective excitation for the coherent state. In fact, the cross 
term, for example, the product operator I^x (/x = x, y, z) (there are also other 
cross terms Ifj,y and /^z in three dimensions), reflects the coupling between the 
atomic center-of-mass and internal motions. In order to excite efficiently the 
coherent state the varying frequency of the force function f{t) should be close 
to the oscillatory frequency lu of the harmonic oscillator. Now one may write 
the unitary propagator J7i(f) for the first-order effective Hamiltonian Hi^{t) of 
Eq. (19) as the form 

UL{t) = exp{-il, [ dt'na{t')}TeM-T [ dt'Hfit')} (20) 
Jo " ^0 

where the first-order effective Hamiltonian Hf{t) is defined by 

Hf{t) = ^ + lmw'x' + fmx (21) 

It is clear from Eq. (3) that the first-order effective Hamiltonian Hf{t) of Eq. 
(21) is very similar to the Hamiltonian (3) of a forced harmonic oscillator. The 
Hamiltonian Hf{t) is responsible for the state-selective excitation of a coherent 
state, as can be seen below. 

It is known that at the initial time of the state-selective trigger pulse [1] the 
halting-qubit atom may be either in the product state |V'o(^)) Iffi) (the atom is in 
the internal state \gi) and the ground motional state |V'o(2;))) or in the product 
state |V'o('^)) Iffo) (or \'tljQ{x))\e)), while the two laser hght beams are applied 
selectively to the two internal states l^o) and |e) and do not have any effect on 
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other atomic internal states including the state |gi). Consider the first case that 
at the initial time the halting-qubit atom is in the product state \ipQ{x))\gi). 
The time evolution process for the atom is expressed as 

\^{x,r,t)) = C/(t)|Vo(x))|gi) = eM-^Eit/n}UL{t)\Mx))\gi). 

Since the atomic operator = (|e)(e| — \go){go\)/2, one has Iz\gi) = 0. Then 
it follows from Eq. (19) that HL{t)\ipo{x))\gi) = iJo|V'o(^)> Iffi)- Note that 
Ho\4>q{x)) = ^tiu)\ipQ{x)) . Thus, the time evolution process is reduced to the 
form 

\^{x,r,t)) = oxp[-tEit/h]e^p[-tujt/2]\^o{x))\gi). 

Up to a global phase factor the state \'if{x,r,t)) is really the initial product 
state |'^Q(3:))|gi). This shows that the Raman laser light beams do not have any 
significant effect on the initial product state \tjjQ{x))\gi). On the other hand, for 
the second case that at the initial time the halting-qubit atom is in the product 
state |'^Q(a;))|g'o) the time evolution process for the atom is expressed as 

|*(a;,r,i)) =exp{-zJ,y^ dt'na{t')}Texp{-^ dt'Hf{t')}\Mx))\9o)- 
Since Iz\go) = "(1/2)150), the wave function is reduced to the form 

|*(x,r,t))=exp{ii^ dt'n,{t')}TeM-'^J^ dt' HJ {t')}\Mx))\9o) 
where the Hamiltonian HJ (t) is given by 

One sees that the Hamiltonian HJ (t) is really the Hamiltonian (3) of a forced 
harmonic oscillator. Thus, the ground state |-0q(x)) will be transferred to a 
coherent state under the Hamiltonian HJ{t) [13, 17]. In an analogous way, one 
can prove that when the halting-qubit atom is in the product state |'(/'g(a;)) |e) 
at the initial time, the time evolution process is given by 

|*(a;,r,t)) =exp{-ii j'^ dt'n,{t')}T e^-'- j'^ dt' H+{t')}\iJo{x))\e) 

where Iz\e) = (l/2)|e) has been used and the Hamiltonian Hj{t) is given by 

Hj{t) = ^ + \muJ-'x^ + \f{t)x. 

The Hamiltonian H^{t) is also the Hamiltonian (3) of a forced harmonic oscil- 
lator and hence the initial ground state |^Q(a;)) will be excited to a coherent 
state under the Hamiltonian Hf{t) [17]. Here gives a summary for the above 
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discussion for the ofF-rcsonancc excitation process. The product state of the 
atom will keep unchanged under the Raman laser light beams when the atom is 
in the product state \tpQ{x))\gi) at the initial time. When the atom is changed 
to the internal state |go) or |e) from the original internal state and hence 
in the product state \tpQ{x))\go) or |^/'o(x))|e) at the starting time of the excita- 
tion process, the atom will be excited to a coherent state by the Raman laser 
light beams and moreover the atomic internal state keeps unchanged during the 
excitation process. Therefore, this coherent-state excitation process is clearly a 
state-selective excitation process. Such a pulse sequence consisting of the two 
Raman laser light beams could act as the state-selective trigger pulse to realize 
the reversible and unitary halting protocol. 

The off-resonance selective excitation using the conventional Raman adia- 
batic laser light beams is a simple technique to prepare approximately a coherent 
state of an atom or atomic ion in the harmonic potential field. The above the- 
oretic investigation for the off-resonance selective excitation shows ones clearly 
the mechanism for a Raman adiabatic pulse sequence to excite selectively the 
coherent state of an atom in the harmonic potential field. The similar theoreti- 
cal analysis also can be seen in a number of references [19, 20, 25a, 31, 33, 45]. 
The state-selective excitation process is indeed involved in the coupling between 
the atomic center-of-mass motion and the internal electronic (or spin) motion 
of the atom, since there is the cross term xl^ in the effective Hamiltonian H]^{t) 
of Eq. (19) which is responsible for the state-selective excitation of the coherent 
state. The coupling between the atomic center-of-mass and internal motions 
is also a general mechanism for the atomic laser cooling [22] and the deceler- 
ating and accelerating processes based on the STIRAP method of a free atom 
[8]. However, if one compares the present theoretic treatment for an atom in 
the harmonic potential well to that one for a free atom [8], one can see their 
difference is quite large. The reason for this is that a motional state of the 
atom in the harmonic potential field is discrete instead of continuous. Though 
the ground motional state of the atom in the harmonic potential well has a 
Gaussian shape, the excited motional state which is generated from the ground 
motional state by the Raman pulse sequence could not have a Gaussian shape 
and tends to have a wave-packet shape quite different from a Gaussian shape. 
For the off-resonance selective excitation above the excited motional state has 
a Gaussian shape in the first-order approximation, but in a general case it does 
not have a Gaussian shape and its shape could be complicated. Thus, in this 
sense the ofF-rcsonance selective excitation mentioned above generally is not an 
optimal technique to construct the state-selective trigger pulse. But it could 
be considered as a starting point to develop further a better technique to gen- 
erate state-selectively the standard coherent state of harmonic oscillator. One 
could exploit the average Hamiltonian theory [37, 39, 54] or the numerical op- 
timization method [38] based on the Heisenberg equations (12) to improve the 
off-resonance excitation technique. The state-selective excitation pulse obtained 
with these optimal methods could be more useful in practice. Of course, it must 
be pointed out that a state-selective excitation pulse is not always equal to a 
state-selective trigger pulse. 
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2.3. The coherent average method 

Here suggests the programming operator composition method to construet 
the state-selective trigger pulse in theory. This method, which has been used 
extensively to construct NMR multi-pulse sequences, is also called the coherent 
average method in nuclear magnetic resonance spectroscopy [54] . The theoreti- 
cal basis of the method is the famous Baker-Campbell-Hausdorff (BCH) formula 
[39, 41] and the Trotter-Suzuki formalism [40, 42]. The BCH formula is also 
the theoretical basis of the average Hamiltonian theory [37, 54]. It is necessary 
to show that the state-selective trigger pulse can be constructed theoretically in 
an error as small as pleased so as to show that both the reversible and unitary 
halting protocol that is insensitive to its input state and the efficient quan- 
tum search process are feasible. Here consider the double-frequency selective 
excitation method, which may be different from the off-resonance selective exci- 
tation method based on the Raman adiabatic laser light beams in the previous 
subsection. The method uses simply two conventional amplitude- and phase- 
modulation laser light beams (the plane wave electromagnetic fields), which 
may not be Raman adiabatic laser light beams, to excite selectively the atomic 
internal states so as to set up the coupling between the center-of-mass and the 
internal motion of the atom. In order to treat conveniently the double- frequency 
state-selective excitation process in the three-level atom system here consider 
only the special case that the two laser light beams have the specific parameter 
settings as given below. It is known the electric dipole interaction between the 
atom and a pair of laser light beams is given by Eq. (8). Now suppose that 
the two laser light beams are amplitude- and phase-modulating such that their 
Rabi frequencies and phases satisfy the following match condition: 

flo{t) = ni{t), v5oW = a + 7, ^i{t) = (wo - iOi)t - a 7. (22) 

Then by using the match condition (22) the electric dipole interaction of Eq. 
(8) is reduced to the form 

Hi{t) = 2hQo{t) exp{—iujot)I^ exp[i^{ko + k\)x — i"f\ cos[i(A:o — k\)x — a] 

+2MlQ{t) exp{iu}ot)I~ exp[— i^(fco -|- ki)x + i'y] cos[^(A;o — ki)x — a] (23) 

where the phases a and 7 can be set suitably in experiment. The theoretical 
treatment could become more convenient for the state-selective excitation pro- 
cess when the electric dipole interaction Hi{t) takes the form of Eq. (23), since 
the electric dipole interaction (23) is modulated by a single frequency uq- On the 
other hand, it could be convenient to treat the state-selective excitation process 
in the rotating reference frame. It is known that the Schrodinger equation for 
a quantum system in the interaction representation can be written as [2] 

in^/i!i{x,t) = Hi{t)^i{x,t) 
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where the state i{x, t) of the quantum system in the interaction representation 
is related to the state (x, in the laboratory frame by the unitary transforma- 
tion: ^i{x, t) = Uo{t)~^^L{x, t), and the Hamiltonian Hi{t) of the quantum sys- 
tem in the interaction representation to the Hamiltonian H{t) = Ho{t) + Hi{t) 
in the laboratory frame by 

Hi{t) = Uo{t)+H,{t)Uo{t) (24) 

where the unitary propagator Uo{t) is defined as 

Uo{t) = Texp{-- dt'Hoit')}. (25) 
The total unitary propagator of the quantum system then is written as 

U{t) = Texp{-^ / dt'[Ho{t')+Hi{t')]} = Uo{t)Ui{t) (26) 
Jo 

where the unitary propagator of the quantum system in the interaction repre- 
sentation is defined as 

Ui{t) = Texp{-| J* dt'Hiit')}. (27) 

One may choose a suitable rotating reference frame or the interaction represen- 
tation for convenient treatment of the state-selective excitation process. First 
of all, the atomic internal Hamiltonian of Eq. (4) is rewritten as H{r) = 
aoE + fujOah + {El — 0-o)\9i) {gi\-, where E is the 3x3 unity matrix. Then 
by neglecting the unity operator the total Hamiltonian (1) for the three-state 
atom system is rewritten as 

H{t) = Ho + fiojah + {El - ao){\gi){gi\) + Hi{t) (28) 

where Hq and Hi{t) are given by Eqs. (2) and (23), respectively. Now the 
rotating reference frame may be defined by the atomic internal Hamiltonian 
Ho{r) = hboolz + {El — aQ){\gi){gi\). In the rotating frame the wave function is 

^i{x,t) = exp[i{Ei - ao){\gi){gi\)t/n]e^{i(JoIzt)^L{x,t), 

and the total Hamiltonian H{t) of Eq. (28) of the atom system is replaced with 
the Hamiltonian Hi{t): 

Hi{t) = — + -mw^x^ + h{uja - ijJo)Iz 
+2Hflo{t)I~^ exp[i-(fco -|- ki)x — i^] cos[-(fco — ki)x — a] 

+2tino{t)I~ exp[-i^{ko + ki)x + i^] cos[i(fco - ki)x - a] (29) 
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where the electric dipole interaction (23) has been used. In the foUowing consider 
a simple case that the on-resonance condition is met, that is, {oja — (jJo) = 0, 
and the amplitude Oo(t) is time-independent, that is, Oo(t) = Oq. Then the 
Hamiltonian Hi{t) of Eq. (29) is reduced to the time- independent form 

Hi{t) = Hi = Ho + Hi{a,^) (30) 

where Hq is still given by Eq. (2) and the electric dipole interaction if/ (a, 7) is 
written as 

Hi{a,'y) = 2fiiloI^ exp[i^(fco + ki)x — 17] cos[^(fco — ki)x — a] 

+2hSloI~ exp[— i^(fco + ki)x + i'j] cos[^(fco — ki)x — a]. (31) 

Obviously, the Hamiltonian Hj{a,j) is dependent upon the phases a and 7. 
One may take the time-independent Hamiltonian Hj of Eq. (30) as the basic 
Hamiltonian to construct the state-selective trigger pulse. The unitary propa- 
gator of the atom system corresponding to the Hamiltonian (30) is written as 

Ui(t) =exp[-iHit/H]. (32) 

On the other hand, in the absence of the two laser light beams the time evolution 
process of the atom system is governed by the Hamiltonian Hq of Eq. (2) and 
its unitary propagator is given by 

Uoit) = exp[-iHot/H], (33a) 

and its inverse propagator by 

U+(t) =eiip[iHot/H]. (33b) 

It can turn out in next section that by using the same Hamiltonian Hq of Eq. 
(2) of the harmonic oscillator one can generate the inverse unitary propagator 
Ug{t) with any time t ^ kw/co up to a global phase factor. In fact, there are 
the unitary operator identities: 

Uo{t)Uo{ti) = Uo{ti)Uo{t) = exp[i/3(ti)]i;, (34) 

where E is the unit operator, exp[2/3(ti)] is a global phase factor, and the time 
ti = 2k'K/u) — t (fc = 1, 2, ...). It follows from Eqs. (33) and (34) that one may 
define the inverse unitary operator U^{t) as 

U+{t) = eM-^Piti)]Uoih). (35) 

Hereafter the unitary operator Uo(ti) is also called the inverse operator of the 
unitary operator Uo{t), although it has a difference of a global phase factor 
from the real inverse unitary operator U^{t). The unitary propagators Uo{t) 
with any time t of the harmonic oscillator in the absence of the two laser light 
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beams can be realized directly in experiment. It follows from Eq. (35) that up 
to a global phase factor the inverse propagator U^{t) also can be realized in 
experiment as the unitary propagator Uo{t\) of the harmonic oscillator can be 
realized in experiment. On the other hand, the unitary propagator C//(t) (32) of 
the harmonic oscillator in the rotating frame also can be realized directly by the 
two laser light beams. These three realizable unitary propagators Uo{t), U^{t), 
and Ui{t) in the rotating frame are the basic unitary propagators to construct 
a general state-selective trigger pulse. 

At the first step the unitary propagator of the electric dipole interaction 
Hj{a,j) is created through the three basic unitary propagators. The unitary 
propagator of the electric dipole interaction Hi{a,"f) is defined by 

[//(a,7,t) — exp[—iHi{a,'^)t/H]. (36) 

It is known from Eq. (30) that the electric dipole interaction Hj{a,j) may be 
written as Hj{a,j) = Hj — Hq. Thus, one may first carry out in experiment a 
composite pulse sequence: 

Uo{ti)Ui{t) = e^p[ip{ti)]U+{t)Ui{t) 

= exp[iP{ti)] exp[iHot/n] exp[-iiJ/i/ft]. (37) 

The physical meaning for the pulse sequence (37) is that a pair of laser light 
beams whose parameters satisfy the match condition (22) and time period is t 
are first applied to the atom in the harmonic potential well, then the pair of 
laser light beams are turned off, and the atom then evolves in the time period 
ti in the harmonic potential well without any external laser light field. If the 
time interval St is sufficiently short, then the unitary operator U^{6t)Ui{6t) 
is approximately equal to the unitary propagator Ui{a,j, St) according to the 
famous Baker-Campbell-Hausdorff (BCH) formula for a product of exponential 
operators [39, 40, 41], 

[//(«, 7, 5t) = U+{St)Ui{St)+0{{6tf). (38) 

Since the composite pulse sequence Uo{ti)Ui{t) in the rotating frame can be 
realized directly in experiment, the composite unitary operator {5t)Ui{5t) 
may be realized up to a global phase factor, as shown in Eq. (37). Thus, the 
unitary propagator Ui{a,j,6t) may be realized in experiment up to a global 
phase factor, as can be seen from Eq. (38). There is an error term 0{{St)'^) 
that is proportional to {St)"^ between the desired unitary propagator Ui{a, 7, St) 
and the composite unitary operator U^{St)Ui{St) on the right-hand side of Eq. 
(38). A much better composition for Uj{a, 7, St) with an arbitrary higher-order 
approximation can also be obtained from the two propagators U^{5t) and U {St) 
by using the Trotter-Suzuki formalism [40, 42]. For example, according to the 
Trotter-Suzuki formalism [42] one can obtain a better symmetric composition 
for the unitary operator Ui{a,-f,St) by the multi-pulse sequence: 

Si{St) = U+{St/2)Ui{St)U+{St/2), (38a) 
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Uiia,j,5t) = SiiSt) + 0{{5tf), 
or a much better symmetric composition by the multi-pulse sequence: 

S2n-l{St) = [S2n-3iPnSt)fS2n-3{{l ' 4p„)5t) [52„-3 (p„<5^)]^ (38b) 

Ui{a,^,5t) = S2n-i{St) + 0{{5tf"+'), 

where n > 2 and p„ = (4 — 4^/^^""^))"^. Thus, without losing generality here 
suppose that the unitary operator Ui{a,j,St) can be constructed as exactly as 
pleased. The unitary operator Ui{a,j, St) will be further used to build up the 
state-selective trigger pulse. 

Now examine the commutation relation between the clcctric-dipole interac- 
tions Hi{a,j) in two different phases 7 = and tt/2. It can turn out from the 
electric dipole interaction of Eq. (31) that the hermite commutation operator 
Q takes the form 

Q = i[Hi{a,0),Hi{a,n/2)] = -IGfi^O^/^ cos^[^(fco - ki)x - a] (39) 

where the commutation relations 2/^ = [/+,/^] and [/+,/+] = [1^,1^] = 
have been used. If the phase a = w/A and the Lamb-Dicke limit is met, that is, 
||(fco — ki)x\\ « 1, then the operator cos^[^(A;o — ki)x — a] can be expanded as 

cos2[i(fco - ki)x -a] = ^ + i(fco - ki)x + O{\\{ko - ki)x\\^). 
Thus, the hermitian operator Q can be written as 

Q = -m^nlh - Sh^{ko - ki)flll^x + 0{\\{ko - ki)x\\^). (40) 

Obviously, the second term on the right-hand side of the operator Q is the cross 
term I^x and is responsible for the state-selective excitation of the coherent 
state. Therefore, at the second step one should construct the unitary propagator 
exp(?A(5). The unitary operator exp(zA(5) can be generated from the unitary 
operators {J7/(a, 7, St)} with the help of the BCH formula [41] and it is realized 
by the multi-pulse sequence: 

Ui{a, 0, St)Ui{a, tt/2, St)Ui{a, 0, St)+Ui{a, 7r/2, St)+ 

= ex.p{-[Hi{a, 0), Hi{a, 7T/2)]{St/n)^} + 0{{St)^). (41a) 

A better composition for the unitary operator exp{iXQ) can also be obtained, 
for example, the following multi-pulse sequence will lead to a better result, 

exp[ig(<5i/?i)2] = Ui{a,0,St/V2)Ui{a,'iT/2,St/V2) 
X [Ui{a, 0, St/V2)+Ui{a, 7r/2, St/V2)+]'^ 

xUi{a,0,St/V2)Ui{a,TT/2,St/V2)+O{{St)'^). (41b) 
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Furthermore, it is possible to obtain a much better composition for the unitary 
propagator exp[iQ{5t/ti)^] by using the Trotter-Suzuki formahsm [42, 43]. Here 
it will not be further discussed in detail. By neglecting the error term 0(||(A;o — 
A;i)x|p) on the right-hand side of Eq. (40) the unitary operator exp[iQ{St/H)'^] 
is written as 

exp[iQ{St/Hf] = exp{-i8fig[4 + AkI^x]{Stf}. (42) 

It can turn out below that the unitary propagator cxTp[iQ{5t/h)'^] can excite di- 
rectly and state-selectively the ground motional state of the harmonic oscillator 
to a Gaussian wave-packet motional state with a high motional energy. 

Now suppose that at the initial time the atom is in the product state cxp(i(/3Q) 
x\ipQ{x))\gi) in the rotating frame, where exp(i(^Q) is a global phase factor. 
Applying the unitary propagator exp[iQ{St/h)'^] to the initial product state one 
obtains 

\^{x,r,t)) = exp[iQ{dt/nf]exp{iipQ)\^o{x))\gi) = exp(iv.o)IV'oW)l5i) (43) 

where the eigen-equation Iz\gi) = 0\9i) has been used. Thus, the initial product 
state keeps unchanged under the imitary propagator cxp[iQ{St/H)'^]. On the 
other hand, if the initial product state of the atom is exp{i(pQ)\'(pQ[x))\go) or 
exp(i(^o)IV'o(^))k) ill the rotating frame, then applying the unitary propagator 
to the initial product state one obtains the atomic product state 

\^{x,r,t)) = exp(iyo) exp{i4nl{Stf} exp{iAD.lAk{dtfx}\ipo{x))\go). (43a) 



or 



|*(x, r, t)) = exp(i(po) exp{-iAnl{6tf} exp{-i'inlAk{6tfx}\ipQ{x))\e). 

(43b) 

where the eigen-equations Iz\9o) = (^l/2)|3o) and Iz\e) = (l/2)|e) have been 
used. If now the rotating frame is changed back to the laboratory frame, then 
only a global phase factor is generated for each of the three atomic product 
states |\I'(x,r, i)) of Eqs. (43), (43a), and (43b). Thus, the atomic product 
states \'^{x,r,t)) of Eqs. (43), (43a), and (43b) are really the final atomic 
product states after the atom is applied by the double- frequency pulse sequence 
that generates the unitary propagator exp[iQ{5t/ti)'^]. Now the ground motional 
state \iPq{x)) of the harmonic oscillator takes the Gaussian form 

where {Ax)q = {2^) and w is the oscillatory frequency of the harmonic os- 
cillator. Obviously, this motional state keeps unchanged during the unitary 
propagator cxp[iQ{6t j K)^] acting on the atom when the atom is in the internal 
state I (71), as shown in Eq. (43). However, the atomic product states (43a) 
and (43b) show that after the unitary operator exp[iQ{5t / hY] acts on the ini- 
tial product state exp(i<^o)IV'o(^))lfi'o) or exp(i<^o)l'^^o(^))k) the atom is in the 
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motional state 



\'i{x,t)) = exp[i^{t)][^--^^-^]^/*exp[--j^-^]exp[ipox/h] (44a) 

or 

\^{x, t)) = exp[i</,(i)][^-^i-^]V4 exp[-i^^] exp[-ipoar/fi] (44b) 

where the global phase factor exp[i(p{t)] = exp(i(^o) exp{i4r2o((5t)^}, exp[i(j){t)] = 
exp(i((9Q) exp{— i4riQ((5t)^}, and the mean momentum po is given by po/H = 
4f2QAA:(5f)^. Both the motional states (44a) and (44b) are the Gaussian wave- 
packet states with the center-of-mass positions a;o = and the momentums 
po and —po, respectively. Thus, the two motional states show that after the 
double-frequency pulse sequence is turned off the atom is still in the original 
harmonic potential well and has approximately the mean motional energy, 

Eq = pl/{2m) + ^nuv'^xl = [AHAknliStf]'^ / {2m) . (45) 

This motional energy Eq can be much larger than the zero-point energy of the 
atom in the harmonic potential well. On the other hand, the motional state 
(44a) shows that at the end of the excitation process of the unitary propagator 
exp[iQ{St/fi)^] the atom in the internal state \go) moves along the direction +x 
with the velocity po/m'm the harmonic potential well, while the motional state 
(44b) shows that the atom in the internal state |e) moves along the direction 
~x with the same velocity po/m. The above investigation for the effect of the 
unitary propagator exp[iQ{St/fi)'^] on the initial product state of the atom is 
summarized as follows. When the atom is in the internal state at the initial 
time, the unitary propagator exp[iQ{6t/fi)'^] does not have any significant effect 
on the initial product state. In particular, even if the initial motional state 
of the atom is an arbitrary wave function, i.e., a superposition motional state, 
the unitary propagator exp[iQ{6t/}i)^] does not yet have any siginificant effect 
on the motional state when the atom is in the internal state \gi) at the initial 
time. However, the initial ground motional state of the atom will be transferred 
to a Gaussian wave-packet state with a higher motional energy by the unitary 
propagator if the atom is in the internal state \go) or |e) at the initial time. 
Therefore, the double-frequency pulse sequences (41) to generate the unitary 
propagator exp[iQ{6t/ti)^] may generally act as the state-selective trigger pulse 
in the quantum control process [1]. 

In order to implement the double-frequency pulse sequences (41) one needs 
to realize not only the unitary propagators {?7/(a,7,t)} but also their in- 
verse propagators {U^{a,'y,t)}. The inverse unitary operators {U/{a,j,5t)} 
can be generated as follows. Since the electric dipolc Hamiltonian Hi{a,j) 
of Eq. (31) is time-independent, these inverse unitary operators may be gen- 
erated by inverting the electric dipole Hamiltonian: Hj{a,'y) — i?/(a,7). 
Since the electric dipole Hamiltonian is dependent upon the phases a and 
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7, it is possible to choose suitably the phase 7 to obtain the negative-sign 
Hamiltonian —Hj{a,j). In fact, it follows from Eq. (31) that Hj{a,0) = 
—Hi{a,Tr) and Hi{a,'jT/2) = —Hi{a, 3Tr/2). Therefore, the inverse unitary op- 
erators U^{a,0,6t) and Uf {a, 1^/2,51) are respectively given by 

Uf{a,0,5t) = Ui{a,'K,8t), U^{a,Tr/2,St) = Ui{a, 3^/2,61). (46) 

Note that the unitary propagator Ui{a, 7, 6t) for any phase values a and 7 can 
be implemented in experiment. Then the unitary propagator exp[i(5(5t/fi)^] 
can be implemented in experiment, since it follows from Eqs. (41b) and (46) 
that the unitary propagator ex])[iQ{St/fi)'^] can be expressed as 

exp[iQ{St/n)^] = Ui{a, 0, dt/V2)Ui{a, n/2, 5t/V2)Ui{a, n, 5t/V2) 

X Ui{a, 37r/2, 6t/V2)Ui{a, tt, St/V2)Ui{a, 37r/2, St/V2) 

xUi{a,0,St/V2)Ui{a,Tr/2,6t/V2) + O{{6t)^). (47) 

The double-frequency pulse sequence (47) of the propagator exp[iQ{dt/fi)'^] con- 
sists of a number of the three realizable imitary propagators Uo{t), U^{t), and 
Ui(t). Obviously, it can be efficiently implemented in experiment, as can be seen 
from Eqs. (37), (38), (38a), and (47). It should be pointed out that one also 
can obtain a much better multi-pulse sequence than the sequence (47) for the 
unitary propagator exp[iQ{St/fi)^] with the help of the Trotter-Suzuki formal- 
ism [42, 43]. This fact tells ones that the unitary propagator exp[iQ{dt/ti)'^] can 
be implemented efficiently and as exactly as pleased. 

On the other hand, the internal-state rotating operator Rz{0) = exp{—i9Iz) 
of the atom in the harmonic potential well in the rotating frame can also be 
prepared in a similar way to generating the unitary propagator exp[iQ((5f/fi)^]. 
When the phase a = and the Lamb-Dicke limit ||(fco — ki)x\ \ « 1 is met, the 
hermitian operator Q of Eq. (39) may be written as 

Q = -i6n^nli, + o{\\iko-h)x\\^). 

If the error term O(||(fco — A;i)a;|p) is neglected, then the unitary propagator 
exp[iQ{St/h)'^] is really the internal-state rotating operator: 

R,{e) = exp[iQ{St/Hf] = exp{-il6nl{Stfh} (48) 

where the rotating angle 6 = 16^1^(61)^. This rotating operator is independent of 
any atomic motional state but applied only to the two atomic internal states l^o) 
and |e) selectively. Here it must be pointed out that the inverse propagator of the 
unitary propagator exp[iQ{6t/h)'^] and the inverse operator of the internal-state 
rotating unitary operator Rz{9) can also be implemented in experiment, as can 
be seen from Eqs. (41) and (47), because both the unitary operator Ui{a,'y,t} 
and its inverse operator {a,j,t) can be implemented in experiment. 

In this section the on-resonance condition = wq has been used in the 
Hamiltonian (29) to simplify the construction of the propagator ex.p[iQ{St/h)'^]. 
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For a general case that the on-resonance condition does not hold, that is, 
i^a 7^ 1^0) one may use tt pulses exp(— iTr/j,) and/or exp(— itt/j,) to refocus the 
term fi{uia — tOo)Iz in the Hamiltonian (29), where the tt pulses exp(— iTr/^) 
and cxp(— ?7r/y) may be generated by an ultrashort laser light pulse. Then in 
the general case the coherent average method can be used as well to build up 
the unitary propagator exTp[iQ{St/fi)'^]. The programming operator composition 
method has also been used to build up the quantum gates and/or the internal- 
state-selective quantum gates in the trapped ion systems [44, 45, 46]. 

3. Manipulating the complex linewidth of a Gaussian wave-packet 
state 

The Gaussian wave-packet motional state or the standard coherent state gen- 
erated by the state-selective trigger pulse must have a higher motional energy 
than the ground state of the harmonic oscillator, but the Gaussian wave-packet 
state generated by the state-selective trigger pulse could not have an expected 
complex linewidth. However, in the quantum control process to realize the 
reversible and unitary state-insensitive halting protocol and the efficient quan- 
tum search process it could be rcqiiired that the Gaussian wave-packet state of 
the halting-qubit atom have an adjustable complex linewidth after the state- 
selective trigger pulse really acts on the atom. Thus, one needs to construct 
a pulse sequence to control the complex linewidth of a Gaussian wave-packet 
state. This pulse sequence combining with the state-selective trigger pulse will 
form a composite state-selective trigger pulse. This composite state-selective 
trigger pulse may manipulate not only the ccnter-of-mass position and momen- 
tum but also the complex linewidth of the Gaussian wave-packet state when the 
halting-qubit atom is really acted on by the composite state-selective trigger 
pulse. Below it is discussed how to generate a pulse sequence to manipulate 
the complex linewidth of a Gaussian wave-packet state. The complex linewidth 
for a Gaussian wave-packet state could be controlled by an external harmonic 
potential field. It is known that the initial motional state of the halting-qubit 
atom is prepared to be the ground state of the harmonic oscillator, which is 
the Gaussian wave-packet state \iI)q{x)) of Eq. (44), and the harmonic potential 
field that is applied to the atom has the oscillatory frequency uj. Now this har- 
monic potential field is switched to another harmonic potential field with the 
oscillatory frequency u)c at the initial time. The atom then undergoes a time 
evolution process in the new harmonic potential field. This evolution process 
may be described by the unitary propagator of the harmonic oscillator with 
the oscillatory frequency Wc, which in the coordinate representation may be 
expressed as [13, 14, 21] 



where the period of the process Tc = h — ta is chosen such that it satisfies 
cos[uJcTc] = and sin[wcT'c] = 1, that is, = [2kTT + Tr/2]/u!c {k = 0, 1, ...). The 




(49) 
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atomic wave function at the end of the process is calculated by 



^{Xb,tb)^ J dXaG{xb,tb;Xa,ta)^oiXa) 

= ^2Mhc)l^'^'\lM Jd-aeM-axl + bx^), 

where the ground state tpQ{xa) of Eq. (44) is used and the parameter a = 
[4(Ax)o]~^ and b = —ih~^muJcXb- By using the Gaussian integral formula: 



/ 





dxa exp(-aa;„ + bxa) = \l i^^) 



one obtains the wave function: 

^{xb,tb) = exp(i'/'o)[^^^(^]'/'exp{-i^}, (51) 

where the phase <j)Q = — 7r/4 and the wave-packet spreading is just e = •\/2(Aa;), 
and (Aa;)^ is given by 

(Ax? = \ - ? = (—)( ^ ) 

^ ' ^2ma;c(Aa;)o ^ujj^2rmjj 

where (Aa;)Q = (5;^) is used. The imaginary part of the complex linewidth of 

the Gaussian wave-packet state (51) is zero and the real part is (A.x)^. There- 
fore, the wave-packet spreading e or the real part of the complex linewidth is 
controlled by both the oscillatory frequencies and w of the harmonic poten- 
tial fields. The initial wave-packet motional state V'o(^) (44) usually may be 
prepared to have a small and fixed wave-packet spreading Eq = •\/2(Ax)o which 
corresponds to a large oscillatory frequency w. For example, if one takes w ~ 10^, 
then (Ax)o ~ 10^^^ and ujh/{2m) ~ 0.1 for the atomic mass m ~ 10~^^Kg. 
Then by setting suitably the oscillatory frequency ujc one may obtain the de- 
sired wave-packet spreading or the real part of the complex linewidth for the 
state ^'(xf,, tb) (51). The important thing is that the center-of-mass position and 
momentum for the Gaussian wave-packet state of the atom keeps unchanged in 
the evolution process, as can be seen from the state Vo(^) (44) and the state 
\E'(xb,t6) (51). This means that after the evolution process the atom is still in 
the original position x = in the coordinate axis and has zero momentum. One 
therefore concludes that the real part of the complex linewidth of a Gaussian 
wave-packet state of an atom may be manipulated by varying the oscillatory 
frequency of the harmonic potential field applying to the atom. 

The imaginary part of the complex linewidth of a Gaussian wave-packet 
state may also be manipulated at will. One of the simplest and most intuitive 
methods to manipulate the imaginary part of the complex linewidth is that the 
atom undergoes simply a free-particle motion or an inverse free-particle motion. 
It is well known that the wave-packet spreading of the Gaussian wave-packet 
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state of a free atom (the initial imaginary part of the complex linewidth is zero) 
becomes larger and larger when the atom undergoes a free-particle motion [2]. 
Obviously, the imaginary part of the complex linewidth will become less and 
less if the atom undergoes the inverse frcc-particlc motion. Suppose^ that at the 
initial time the atom in the harmonic potential well is in the Gaussian wave- 
packet state of Eq. (51). Now one turns off the harmonic potential field applying 
to the atom. Then the atom becomes a free atom. However, the atom does not 
leave its original position even after the harmonic potential field is switched 
off, since the atomic motional momentum is zero before the harmonic potential 
field is turned off. Therefore, the only change for the Gaussian wave-packet 
state (51) of the free atom is its complex linewidth after the harmonic potential 
field is turned off. The time evolution process of the atom after the harmonic 
potential field is turned off may be calculated by using the free-particle unitary 
propagator. It is known that the unitary propagator of a free particle is given 
by [13, 14, 21] 



G{xb,U;Xa,ta) = y:j^:^^Mi^ixb-Xaf] (52a) 

where T = tf, — ta is the time period of the free-particle motion. After the 
harmonic potential field is turned off the free atom undergoes a free-particle 
motion with the time period T and its motional state 'if{xa,ta) of Eq. (51) is 
changed to another Gaussian wave-packet state: 

Here the complex linewidth of the Gaussian wave-packet state f{xb,tb) is 
given by W{T) = (Aa;)^ -|- iHT/{2m). Thus, the imaginary part of the com- 
plex linewidth is proportional to the time period T of the free-particle motion. 
However, the imaginary part hT/{2m) is always positive. In order to achieve a 
negative imaginary part one may let the atom perform an inverse free-particle 
motion. The unitary propagator for the inverse free-particle motion may be 
given by 

G+{xb,tb;Xa,ta) = J - .^^^^ exp [- i — [Xb-Xaf] (52b) 

The unitary propagator of the inverse free-particle motion may be generated 
with the help of the external quadratic potential field (see below). Now the state 
^{xa,ta) of Eq. (51) is changed to the Gaussian wave-packet state ^i{xb,tb) 
after the inverse free-particle motion, 

^.fa. = explVol l^]-/' exp{-i ^^^^/j ^^g^ }. (53b) 

Here the complex linewidth of the state 'i>i(xb,tb) is given hy W = (Ax)^ — 
ihT/{2m). Its imaginary part is negative. Thus, the imaginary part of the 
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complex lincwidth can be controlled by the free-particle motion and the in- 
verse free-particle motion. Obviously, due to the fact that the atomic motional 
momentum is zero the center-of-mass position {xq = 0) of the Gaussian wave- 
packet state (51) keeps unchanged in the frce-particlc motion and the inverse 
free-particle motion. The manipulating method by using the free-particle mo- 
tion and the inverse free-particle motion is very simple. However, there could be 
a disadvantage for the manipulation: if a large imaginary part of the complex 
linewidth needs to be achieved, then one needs to spend a long time for it. 

In the preceding discussion one needs to use the inverse unitary propagator 
of a free particle to manipulate the imaginary part of the complex lincwidth. 
Here gives the explicit expression for the inverse unitary propagator without 
a detail proof. It is known that the unitary propagator of a free particle is 
given by Uj{t) = cxp[— ^^], p is the momentum operator of the free particle. 
Denote that Uok{tk) = exp[— iifofe^fe/^] = 1)2) is a unitary propagator of a 
harmonic oscillator. The Hamiltonian iJofe of the harmonic oscillator is i/ofe = 
/{2m) + muf.x'^ /2. Then it can prove that the inverse imitary propagator of 
a free particle can be written as, up to a global phase factor, 

Uf{T)+ ^ exp[z^] = Uo,{T,)UfiT)Uo2iT2) 

where the time intervals Ti and T2 (or the oscillatory frequencies ui and 0J2 of 
the two harmonic potential fields) are determined through 

sm(a;iri) = =F- 



and 





a;^)]2 + [2Ta;ia;i]2 ' 








- a;2)]2 + [2Ta;iwi]2 


2Tuj 




^\T^u:\uj\ + (^i - 


c^2)]2 + \2Tu:\u:2f ' 







cos(wiTi) = =F 
sin(a;2T'2) = ± 

cos(a;2T2) = ± , „_ „ „ „ „ 

^ ^ ^/{T''u:lJi + (a;i - a;?)]2 + \2Tw\u-,\^ 

Given the time period T of the free-particle motion and the oscillatory frequen- 
cies LO\ and L01 of the two harmonic potential fields one can calculate the time 
intervals T\ and T2 from the four equations above. One sees that the inverse 
unitary propagator of a free particle can be realized only when the specific 
external harmonic potential fields are applied to the particle. Therefore, the 
time evolution process of a free particle is a self-irreversible evolution process, 
although this process is unitary. 

To improve the simple method based on the free-particle motion or the in- 
verse free-particle motion more complex multi-pulse sequences may be employed 
to adjust the complex linewidth. The pulse sequences consist of several pulses 
of the harmonic potential fields with different oscillatory frequencies. One of the 
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pulse sequences is given below. It is known that after the free-particle motion 
with a short time period T the atom is in the Gaussian wave-packet state of Eq. 
(53a). Now the atom is applied by a pulse sequence consisting of two harmonic 
potential field pulses with different oscillatory frequencies (see below) such that 
the time evolution process of the atom is described by the unitary propagator; 



G{xb,tb,Xa,ta) = y^^^^^|^;j^^exp{i(^)[S'66a;j -hS'a62a;aa;6]} (54) 

where these parameters in the propagator are given later. Now one has the 
initial state of Eq. (53a) with the time parameter T = Tq and the propagator of 
Eq. (54). Then one can determine the time evolution process of the atom under 
the harmonic-potential-field pulse sequence. The final state of the process is 
given by 



'^{xb,tb) = exp(i(/)o)[ 



27r V iH 



X exp{-i(^)[(^!^^)5„\ + ^inSl, - Sbb)]xl}. (55) 
The Gaussian wave-packet state \E'(a;b,f(,) (55) has the complex linewidth: 

W = ( ^ ) ^"^"^ ~ ~ (56) 

where S = 2m,{Ax)-^/H. Denote the complex linewidth W = {Ay)"^ + iHTyj/{2m). 
Then one has 

(TpS^f, - Sbb) 
6^St, + {ToSl,-Sbby 



and 

^'Si, + {ToSj,-Sbb)^ 



i^yf = .... ^^f^' . (58) 



Obviously, > if {ToSlt, - Sbb) < and < if {ToSl, - Sbb) > 0. Only 
when S'^S^f^ + (ToS^f, — Sbb)'^ «1 can it be possible for the absolute parameter 
\Tw\ to be much larger than one. The equations (57) and (58) show that both 
the real ((Ay)^) and imaginary (T^) parts of the complex linewidth can be 
controlled simultaneously by the harmonic-potential-field pulse sequence. 

Suppose that the harmonic-potential-field pulse sequence is applied to the 
atom in the manner that the first harmonic potential field with the oscillatory 
frequency u) is applied to the atom at the initial time, it lasts a time interval T, 
then it is turned off and at the same time the second harmonic potential field 
with the oscillatory frequency lOq is turned on, and then it lasts a time interval 
To. Thus, the total time period of the harmonic-potential- field pulse sequence 
is T + To- It is known that the unitary propagator of a harmonic oscillator is 
generally written as [13, 14, 21], 



G(xb,t,;x„,M = J ,2^^,i,(^r) ^"P^^ 2?isin(a;T) ^^^^ + ^"^^^"^^ " '^^^"^^ 

(59) 
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The propagator (49) is a special case of the propagator (59) with the time 
period T = [Ikn + 7r/2]/a;. One can see that when the time period T satisfies 
ujT = kn, the propagator (59) appears singular as sin(ci;T) = 0. It seems that 
the propagator of the harmonic oscillator should take a different form from the 
original one (59) at these time period points that make the propagator singular. 
But it can turn out that even at these time period points the correct propagator 
of the harmonic oscillator can be obtained directly from the original one (59) 
[47, 48]. Therefore, the unitary propagator (59) indeed can describe completely 
the time evolution process of a harmonic oscillator with any time period. The 
composite unitary propagator of the harmonic-potential-field pulse sequence 
then is calculated by 

G{Xb,tb;Xa,ta) = j dXcG2{Xb,tb;Xc,tc)Gi{Xc,tc;Xa,ta) (60) 

where {Gk{x' ,t';x,t), k = 1, 2} are the unitary propagators of the atom under 
the first {k = 1) and the second (k = 2) harmonic potential field, respectively. 
Note that T = tc - ta, To = h - tc, and U - ta = T + To. By substituting the 
unitary propagator (59) of the harmonic oscillator in the equation (60) and then 
by a complex calculation one obtains 



G{Xb,tb,Xa,ta) = ^ ( ) exp{i{^)[Sbbxi + Sab^XaXb + SaaXl]} (61) 

where the parameters are given by 

rj = [ujo cos(uJoTo) sin(cL'T) + Losi'n{LOoTo) cos(a;r)], 

Sbb = —[—uJoSm(u!T)sin(LL)oTo) + to cosiuoTo) cos(ujT)], 
V 

'~>ab — ) 
V 

Saa = — \—ujsm.{ujT)sm.UjOoTo) + uJ o cosiujT) cosiijj oTg)]. 
V 

If one sets the parameter Saa in Eq. (61), then the propagator (61) is 
reduced to the propagator (54). Since the parameter r] satisfies < |ry| < ui+cOo, 
the parameter Saa = means that 

tan(a;r)tan(a;oro) = — . (62) 

UJ 

By using the equation (62) one can reduce respectively the parameters Sbb and 
S^i^ to the forms: 

_ {u^ - tol) tan(a>r) 

^bb - ■ r, , , _ 9/ 



s: 



2 



oj [l+ta.n^{uiT)]' 

UJo + iMl'l^'T) 



l + tan2(a;r) 
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From these parameters one sees that both the wave-packet spreading \/2(Aj/)^ 
(58) and the time interval T^, (57) are dependent upon the parameters w, T, 
ujo, and To of the two harmonic potential fields. There are four undetermined 
parameters w, T, uJq. and To of the two harmonic potential fields, while there 
are only three independent equations (57), (58), and (62) to determine these 
parameters. Thus, given the complex linewidth W = {Ay)'^+iHTi^/ (2m) for the 
Gaussian wave-packet state ^'(.Xb.tb) (55), one can determine these parameters 
for the two harmonic potential fields. 

The equations (57) and (58) lead to the relation: 

(Ax)^T^ _ {TpSl, - Sj,,) 

{^y? ~ SI, ■ ^^'^ 

Substituting the parameters and 5'^^ in the equation (63) one obtains 

ujT^{/\xf _ {u:TQ)[nl+t&n^{ujT)\ - (1 - nl)t&n{LoT) 



[Ay? [n2+tan2(wT)] 



(64) 



where = oj1/lo^. Now the equation (64) and the parameters S},b and S'^;, are 
used to simplify Eq. (58) to the form 

u^^Axf 2m{Ayf [1 + tan^(a;T)] 

lA^^^ + " [n2+tan2(a;T)]- ^^^^ 

Therefore, one obtains these three independent equations (62), (64), and (65) 
which can be used to determine the three independent parameters rio, tOoTg, and 
ujT if one is given in advance the parameters (Ay)^, (Aa;)^, Tyj, co, and coTq. 
First of all, one can solve Eq. (65) to obtain the parameter : 

, ^ 1 + {1 - (^Tg^)[(^^)^ + T^]} tan^(^T) 

^^o . t^2(Aa;)2 2m(Ay)^ ^2 i ^21 ' ^ ' 

I (Ay)^ )i\ n +^w\ 

Then inserting the parameter into Eq. (64) one obtains 

tan(a;T) = - ^ (67) 

where the parameters Aq and Bq are obtained from the given parameters (Ay)^, 
(Aa;)^, Tyj, u, and wTq through the equations: 

^~(Ayf~ — ft — 



5o=a;To + ^^^±^^. (68b) 



wT^(Aa;)2 
"(A^ 

Now the parameter ojT can be conveniently determined from Eq. (67) when 
the parameters Aq and Bq are obtained in advance. After the parameter wT 
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is obtained one can determine the parameter Hq = ujo/<-^ from Eq. (66) and 
further obtain the parameter Wo^o from Eq. (62) by using the parameters u>T 
and Uo- Obviously, given different parameter values (Ay)'^ and T^, while the 
other parameters (Aa;)^, uj, and ojTq arc kept constant, one can obtain a different 
parameter set {ojo, To, uj, T}. The set of parameters then are used to generate 
the two harmonic potential field pulses. 

As a summary, in the above discussion one first uses a harmonic poten- 
tial field to adjust the parameter (Aa;)^ and then uses a free-particle motional 
process and a pair of harmonic potential field pulses to adjust jointly the com- 
plex linewidth of a Gaussian wave-packet motional state. In these processes 
the center-of-mass position (xq = 0) and momentum (po = 0) of the atomic 
Gaussian wave-packet motional state always keep unchanged due to the fact 
that the atomic motional momentum is zero. Finally, it can turn out that the 
inverse unitary propagator of a harmonic oscillator may be prepared by its uni- 
tary propagator. It is well known that the inverse propagator U{t, to)"*" of a 
unitary propagator U (t, to) is just equal to the unitary propagator U {Iq, t), that 
is, U(to,t) = U{t,to)~^. It is known that the unitary propagator of a harmonic 
oscillator is given by Eq. (59). If one sets the time interval of the propagator 
(59) toheT = th-ta = 2kTv/uj - T' or uT = 2k-K - uT', then the propagator 
(59) is rewritten as 



By comparing Eq. (69) to Eq. (59) one sees that up to a global phase factor 

the unitary propagator U{T) (69) is really the inverse propagator U (T")+ of the 
harmonic oscillator with the time period T', that is, U{T) = exp(i</)Q)[/(T')+ 
with a global phase factor exp(i^o)- This means that the inverse unitary prop- 
agator of a harmonic oscillator can be generated from its unitary propagator 
(59). One therefore concludes that the Hamiltonian Hq = /{2m) + muj'^x^ /2 
of a harmonic oscillator can generate both the unitary propagator U{T') and its 
inverse propagator U{T') of the harmonic oscillator. This is completely differ- 
ent from the case of a free particle. As shown in Eq. (38) and (38a) , one needs 
to use the inverse unitary propagator of a harmonic oscillator to build up the 
state-selective trigger pulse. 

4. Manipulating a Gaussian wave-packet state by the unitary prop- 
agator of a general quadratic Hamiltonian 

In the section it is investigated in detail how the unitary propagator of a 
quadratic Hamiltonian (or Lagrangian) affects an atomic Gaussian wave-packet 
motional state in an internal-statc-independent form. A general method to 
manipulate a Gaussian wave-packet state of an atom is to use the unitary prop- 
agator generated by a quadratic Hamiltonian of the atom. A unitary propagator 
generated by a quadratic Hamiltonian does not change the Gaussian shape of 




X exp{— i 



2nsm{wT') 



[{xl+xl)cos{ujT')-2xbXa]]. 



(69) 
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a Gaussian wave-packet state to any other shape when it acts on the Gaus- 
sian wave-packet state [5, 6, 13, 14]. Generally, a quadratic Hamiltonian of a 
quantum system in one dimension may be written as 

H{t) = ^p' + V{x,t). (70) 

Here the generalized quadratic potential operator V{x, t) consists of only linear 
and quadratic terms of the center-of-mass coordinate and momentum of the 
quantum system, 

V{x, t) = h){t){px + xp) + ic(i)a;2 + d{t)p + f{t)x. (71) 

Several typical examples have been given in the previous sections. It is well 
known that the potential operator V{x,t) = for a free particle, V{x,t) = 
muj{tYx'^ /2 for a harmonic oscillator, and V{x,t) = muj{t)'^x^ /2 + f{t)x for a 
forced harmonic oscillator. In the potential operator (71) the linear terms are 
only responsible for manipulating the center-of-mass position and momentum of 
a Gaussian wave-packet state, while the quadratic terms can be used to control 
the complex linewidth of a Gaussian wave-packet state. It is well known that in 
the coordinate representation the unitary propagator of a quadratic Hamiltonian 
(or Lagrangian) can be exactly obtained by the Feynman path integration [13, 
14, 21]. The time evolution behavior of a quantum system with a quadratic 
Hamiltonian has been studied extensively and thoroughly [5, 6, 10, 11, 13, 14, 
21, 47, 48, 49]. The unitary propagator of a quantum system with a quadratic 
Hamiltonian in one- dimensional coordinate space may be generally written as 
[13, 14, 21, 49] 



/ Tfl 

G{Xb,tb;Xa,ta) = \ exp{i^[Sbbxf, + Sab^XaXb + SaaXi]} 

V ilmjab 

i i 
X exp{+-XaQa{tb, to) + jXbQbiU, to)} cxp[je(4, to)] (72) 

where the function fab = {—Sab)~^ [49]. Some important and frequently using 
unitary propagators which are the special forms of Eq. (72) have been given 
in the previous sections: (?') the Hamiltonian of a free particle is H = j [2m) 
and the unitary propagator is given by Eq. (52a); {ii) the Hamiltonian of a 
harmonic oscillator is given by Hq of Eq. (2) and the unitary propagator is 
given by Eq. (59): {Hi) a forced harmonic oscillator has the Hamiltonian of Eq. 
(3) and its unitary propagator is given by [13, 14, 21] 

Gf{xb, h; Xa, to) = exp{iQ{tb, ta)}G{xb, h; Xa,ta) 



X exp{-[Qb{tb,ta)Xb + Qa{tb,ta)Xa]} (73) 
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where G{xb, h; Xa, ta) is given by Eq. (59) with the time period T = tb — ta and 




1 



f(t) sm{u:(tb - t)]dt. 



sin(a;T) 



t. 




-L-jyit)sinHt-ta)]dt, 



muiH sm.{ijjT) 



1 




f{t)f{s) sin[u}{tb — t)] sin[u;(s — ta)]dsdt, 



It has been shown that the unitary propagator (73) of a forced harmonic oscil- 
lator can be used to generate a standard coherent state of a harmonic oscillator 
[13, 14, 15, 16, 17]. It is known in the section 2 that the harmonic potential 
field {muj{t)^x'^ /2) and the forced field {f{t)x) for a harmonic oscillator can be 
generated by the external driving electric or magnetic field. 

If one names the unitary propagator (72) of a quadratic Hamiltonian the 
quadratic unitary propagator, then it can turn out that a product of any two 
quadratic unitary propagators is still a quadratic unitary propagator. This 
property is particularly important as it leads to that a complex quadratic uni- 
tary propagator may be decomposed into a sequence of simple quadratic unitary 
propagators. This makes it convenient to implement a complex quadratic uni- 
tary propagator in experiment, since a simple quadratic unitary propagator can 
be prepared easily in experiment. The unitary propagator (72) consists of the 
quadratic terms x^, 2xaXb, and x^ and the linear terms Xb and Xa in addition to 
the global phase Q{tb, ta)- Suppose that there are two quadratic unitary propa- 
gators, each of which contains only quadratic terms x^, 2xaXb, and x^. Then it 
can prove that a product of the two unitary propagators is still a quadratic uni- 
tary propagator that contains only the quadratic terms. This is a direct result 
of the Lie group generated by the Lie algebra su{l, 1) whose three basis ele- 
ments may be taken as p^, x^, and {px + xp)/2. On the other hand, a quadratic 
unitary propagator which contains linear terms times another quadratic unitary 
propagator that contains only quadratic terms will generate a quadratic unitary 
propagator that contains linear terms. The quadratic terms of the unitary prop- 
agator (72) can manipulate the complex lincwidth of a Gaussian wave-packet 
state, while the linear terms are used to control only the center-of-mass position 
and momentum of the Gaussian wave-packet state. This will be proven below. 

It is known that a standard Gaussian wave-packet state of an atom with 
mass m may be written as [2, 3, 4, 5, 6, 7, 8] 



*o(a;,i) = exp{i(j)Q)[ 



i^x)\ l/4 



1 



2^J 



[(Ax)2+i(|^)] 



xexp{-- 



4[{Ax)^+i{!^)] 



} ex.p[-ipox/fi] 



(74) 
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where xq and —po are the center-of-mass position and momentum of the Gaus- 
sian wave-packet state, respectively, exp(i(/)o) is a global phase factor, and the 
complex linewidth of the Gaussian wave-packet state is given by 



2 , ,/fi^O^ 



W =(Axy + i(-^). 

The Gaussian wave-packet state ^o{x,t) (74) has the wave-packet spreading 
e{To) = Y^2[(Ax)M-~(^^^^y)^. The relation between the wave-packet spread- 
ing and the complex linewidth is given by 

\Wf = ^{Axfsinf. (75) 

The physical meaning for the Gaussian wave-packet state '^o{x,t) is clear: if a 
free atom is in the Gaussian wave-packet state, then the Gaussian wave-packet 
state tells ones that the atom moves along the direction x in the coordinate 
axis with the motional velocity {—po/m). A Gaussian wave-packet state is com- 
pletely described by the three parameters: the center-of-mass position xo, the 
mean momentum (— po), and the complex linewidth W. Now examine the time 
evolution process of the Gaussian wave-packet state ^'o(a;a,<a) (74) under the 
action of the unitary propagator G{xb,tb;Xa,ta) (72) of a quadratic Hamilto- 
nian. The time evolution process may be calculated by 

'^{Xb,tb) = j dXaG{Xb,tb;Xa,ta)^oiXa,ta) 





/ m 


[(Ax)2+i(|^)]\ 


' i^nhfab 



X / dXa{exp{iSc/n} exp{- ^ ^"''^ "^"^ } exp{-ipoXa/n}} (76a) 
J-oo ^[[ixx) +i[2;^)\ 

where the action Sc can be found from the propagator (72), 

5c = ^m[Sbbxl + Sab2XaXb + SaaXl] + X^Qaih, to) + XbQbih, to) + ^©(^6, ^a)- 

One can write the final state ^{xb, h) as, with the help of the Gaussian integral 
(50), 





1 m 


[{Axy + i{^)]\ 


1 i^-Khfab 



(76b) 
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where the parameters a and b are given by 

TTl 1 1 

and 

6 = -S.,x, + -[g„(i„t„) - po] + 2 [(A.)^+,(^)] - 

Since 6^/(4a) is quadratic with the coordinate Xb, the wave function '^{xb,t},) 
of Eq. (76b) can be written as 

^{xb,h) = exp(i(/)o)exp[ie(t6,i«)][^^^]i/^ 



y /a6*aa(Ax)^ + l( 3^) [1 + To^aaJ 

where the complex coefficients A, B, and C are determined through 

+ s.. + c = ^ + '^Sbbxl + '-xMtM - ip^^^Tlf] ■ ^''^ 

By substituting the parameter a of Eq. (77a) and h of Eq. (77b) in the equation 
(79) one can obtain exphcitly the complex coefficients A, B, and C. Furthermore 
one can find, by a complex calculation, 

exp{Re(A)xg + Rc(S).Tb + Rc(C)} 

Therefore, the wave function '^{xbjh) of Eq. (78) can be further written as 
^{xb, tb) = exp(i^o) exp[ie(t6, to)] exp{ j Im(C)} 



X 



-(A^,i/4 / [-Sui>) 



X exp{i[Im(A)a;g + lm.{B)xb]}. (80) 

The wave function ^{xb, tb) of Eq. (80) is indeed a Gaussian wave-packet state. 
This can be seen more clearly from its probability density \^{xb,tb)\'^ which is 
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a standard Gaussian function. Thus, the Gaussian wave-packet state ^{xb,tb) 
(80) has the center-of-mass position: 

= -S^t, i^aaXo + ( )] 

and the wave-packet spreading: 



The imaginary parts of the coefficients A, B, and G in the state 'i{xb,tb) (80) 
are given by 

T^/ A\ _ /'^^ 'S'6b(l + SggTo) + Z{To){SbbSaa - 5^fc) 

Im(B) = c"'^-) (^fc + SabXo){l + SaaTn) + Z{Tn){VbSaa - VgSab) , , 

^ ^fi' [l + SaaToV+5'S^gg ' ^ ^ 

lm(C) = ( ) ^^O*^"" ~ ^O^a)[l + SggTo] - S'^SggVl 

' ^2fi^ (^2^2^ + [1 + 5„„ro]2 ' ^ 

where Z{To) = [SggS^ + Tq{1 + SggTo)] and the parameters S, Vg, and Vb are 
defined by 

2m(Aa;)^ 1 f 

S={ ^ ), Vg= —[Qg{tb,tg) - PO], Vb = —Q b{tb ,tg) . 

The state ^{xb,tb) (80) can be further written as a standard Gaussian wave- 
packet state: 



*(x,i) = exp(z^,)[i^]V4 / 1 



27r Y[(Ay)2 + i(^)] 



' "^^"^ [(Ay)-+t4)] ^"^'"''^'^''^ ^''^ 
where exp(i(^Q) is a global phase factor, the mean momentum (— p) is given 
through 

p/n = - Im(S) - 2xc Im(A), (83) 
the real part of the complex lincwidth is 

, ^ Slb{Slg{^xf + {^ni + SggT,Y} 

Sl + l&{\ra{A)Y{Slg{Axf + {^^ni + SggT,fY' > 
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and the imaginary part equals 



( 



nx, 



w 



) 



st,+mm{A)r{suAxr + { 



Um.{A){Sl,iXrf + ( 



)-[i + ^s'„„ro]-"} 



(85) 



2m 



)2[l + 5„„To]2}2 



2m{Ax) 



Here the imaginary part of the complex linewidth can be positive or negative, 
which is dependent on the term Im(A). It is known from Eq. (81a) that the term 
Im{A) is dependent on only the parameters Sbi,, Sab, and Saa of the quadratic 
terms in the unitary propagator (72) of the quadratic Hamiltonian (70) but in- 
dependent of those parameters of the linear terms. Then the real ((Ay)^) and 
the imaginary (T^) part of the complex linewidth are dependent on only the 
quadratic terms but not the linear terms of the unitary propagator (72). One 
therefore conclndcs that the complex linewidth of a Gaussian wave-packet state 
can be adjusted only by the quadratic terms in the unitary propagator (72). 
This is a convenient method to manipulate the complex linewidth of a Gaussian 
wave-packet state by adjusting only the quadratic terms of the unitary propa- 
gator (72). In the quadratic Hamiltonian (70) the linear terms d{t)p and f{t)x 
do not have a contribution to the quadratic terms of the unitary propagator 
(72) [13, 14, 21, 49]. Then one can use only the quadratic operator terms of the 
Hamiltonian (70) to adjust the quadratic terms of the unitary propagator (72). 
Consequently one can manipulate at will the complex linewidth of a Gaussian 
wave-packet state by the quadratic operator terms of the quadratic Hamiltonian 
(70). On the other hand, it is known from Eq. (81b) that the term lm{B) is 
dependent on the parameters Qaih, ta) and Qbih, ta) of the linear terms of the 
unitary propagator (72). Then the center-of-mass position Xc and momentum 
{—p) of the Gaussian wave-packet state (82) are dependent on the parameters 
Qa{tb,ta) and/or Qb{tb,ta), although they are also dependent on those param- 
eters of the quadratic terms. One therefore can manipulate the center-of-mass 
position and momentum of a Gaussian wave-packet state through the linear 
terms of the unitary propagator (72) or more conveniently through the linear 
terms of the quadratic Hamiltonian (70). A typical example is that in order to 
generate a standard coherent state of harmonic oscillator with a higher motional 
energy one may use the linear term f{t)x [13, 14, 15, 16, 17, 21, 49] that may be 
generated by the external driving field, while one may use the time-dependent 
and frequency-varying harmonic potential field to adjust the complex linewidth 
of a Gaussian wave-packet state, as shown in the previous section 3. 

5. Discussion 

The state-selective trigger pulse has the two basic properties. One of which is 
that the state-selective trigger pulse can have a real effect on the atom only when 
the atom is in some given internal states. Another is that the state-selective trig- 
ger pulse does not change Gaussian shape of an atomic Gaussian wave-packet 
motional state to any other shape. The former property is inherent and the last 
one is due to the fact that a Gaussian wave-packet motional state for a single 
atom is simple and easy to be manipulated and controlled in time and space. 
An internal-state-dependent selective excitation process of an atomic system is 
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generally closely related to manipulation and control in time and space of the 
atomic internal electronic (or spin) motion, the atomic center-of-mass motion, 
and the coupling between the internal and the center-of-mass motion, although 
an inhomogcncous external magnetic field could also generate an intcrnal-state- 
dependent force exerted on a spin and hence could be used to generate in an 
internal-state selective form a coherent state of the spin in the harmonic poten- 
tial field. Therefore, a general construction for the state-selective trigger pulse 
is generally involved in using the electromagnetic field pulses, i.e., the laser light 
pulses, to create the interaction between the atomic internal states and center- 
of-mass motional states and realize the coupling between the centcr-of-mass and 
the internal motion of the atom. A state-selective trigger pulse transfers one 
atomic Gaussian wave-packet motional state to another with the help of the 
atomic internal states. This point is the same as those of the unitary deceler- 
ating and accelerating processes [8]. On the other hand, there are many other 
space-dependent processes in quantum information science. Typical examples 
include the quantum communication [52] and the construction of quantum gates 
by the short-range interactions such as the conditional collision interaction [53] . 
These space-dependent processes generally use the particle (photon or atom) 
transport process to realize the quantum state transfer of the atomic internal 
states or the photon polarization states in space and implement the quantum 
gate operations of the internal-state quantum bits. Here motional states are 
usually considered as carrier of the quantum information transport. These pro- 
cesses generally emphasize realization of the internal-state transfer in space (i.e., 
quantum information transfer) or the quantum gate operations of the atomic 
internal-state quantum bits with the help of the motional states of the parti- 
cles instead of the unitary manipulation of the motional states themselves. In 
this sense these space-dependent processes are different from the state-selective 
excitation process of the trigger pulse and the unitary decelerating and acceler- 
ating processes. In the reversible and unitary halting protocol [1] the quantum 
program converts the difference of the atomic internal states (the initial func- 
tional states) into the difference between the wave-packet motional states of the 
halting-qubit atom in space. Thus, one has to manipulate in time and space 
the atomic wave-packet motional states, in order that the reversible and unitary 
halting protocol is state-insensitive. Note that here the atomic wave-packet mo- 
tional states are not used as quantum bits and the atomic internal states still 
act as the halting quantum bit. 

Manipulating unitarily in space a wave-packet motional state of an atom or 
a superposition of motional states of the atom is generally difficult in exper- 
iment with respect to manipulating a purely time-dependent quantum state. 
However, it is of crucial importance to realize both the reversible and unitary 
state-insensitive halting protocol and the efficient quantum search process. The 
quantum control process that simulates the reversible and unitary halting proto- 
col contains the conventional particle transport process such as the free-particle 
motion, but the more important is that it also contains the coherent-state se- 
lective excitation process of the trigger pulse and the unitary decelerating and 
accelerating processes that are different from the conventional particle transport 
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process. The conventional transport process usually need not require the coop- 
eration of the center-of-mass and the internal motion of the atom, while these 
unitary processes need to manipulate not only the atomic internal and center- 
of-mass motions in space but also their coupling. Unitary manipulation for the 
atomic internal states which act as quantum bits is an important research area 
in quantum computation in the past decade, and it is usually more convenient 
than for atomic wave-packet motional states. One of the main reasons why it is 
generally difficult to manipulate at will an atomic wave-packet motional state in 
space is that the unitary dynamical process is generally complicated for an atom 
in a general potential field. In few simplest potential fields such as a harmonic 
potential field the quantum dynamical behavior of an atom can be completely 
understood, while a complete knowledge for an atom system and its quantum 
dynamics is the basis to manipulate at will the atomic motional states in space. 
Thus, at present the unitary manipulation in time and space of motional states 
of the halting-qubit atom in the quantum control process has to be restricted 
to a simplest case such as the atom in a Gaussian wave-packet motional state 
and in a harmonic potential field. In contrast to the preparation of quantum 
gate operations, here the unitary manipulation in time and space of the atomic 
motional states becomes the main research area. 

The quantum dynamics of a harmonic oscillator and a Gaussian wave-packet 
state have been studied extensively and thoroughly in quantum mechanics. A 
Gaussian wave-packet state is one of the simplest quantum wave-packet states 
that can be manipulated and controlled in time and space easily and precisely. 
In the previous paper [8] it has been shown that the unitary decelerating and 
accelerating processes based on the STIRAP method can transfer one Gaus- 
sian wave-packet motional state of a free atom to another in the ideal or near 
ideal adiabatic condition and can manipulate the center-of-mass position and 
momentum of a Gaussian wave-packet motional state. The advantage of the 
manipulation is that the space-selective manipulation can be carried out easily, 
since the manipulation uses the STIRAP laser light pulse sequence, while laser 
light makes it easy to perform the space-selective and internal-state-selective 
operations of an atom. However, the shortcoming of the manipulation is that 
the complex linewidth of a Gaussian wave-packet motional state can not be 
manipulated at will by the STIRAP method. Now in this paper several meth- 
ods have been developed to manipulate the complex linewidth of a Gaussian 
wave-packet motional state of an atom. Their basic starting point is to apply a 
quadratic potential field to the atom. These results in the previous [8] and the 
present paper show that a Gaussian wave-packet motional state of an atom can 
be manipulated at will in experiment. It can be predicted from these results 
that there is no longer unsurpassable obstacle in theory for a quantum computer 
to solve the unsorted quantum search problem in polynomial time. 
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